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APPLICATIONS OF MODERN HYDRODYNAMICS TO AERONAUTICS. 


By L. PRANDTL. 


PART I. 
FUNDAMENTAL CONCEPTS AND THE MOST IMPORTANT THEOREMS. 


1. All actual fluids show internal friction (viscosity), yet the fofces due to viscosity, with 
the dimensions and velocities ordinarily occurring in practice, are so very small in comparison 
with the forces due to inertia, for water as well as for air, that we seem justified, as a first ap- 
proximation, in entirely neglecting viscosity. Since the consideration of viscosity in the 
mathematical treatment of the problem introduces difficulties which have so far been overcome 
only in a few specially simple cases, we are forced to neglect entirely internal friction unless we 
wish to do without the mathematical treatment. f l 

We must now ask how far this is allowable for actual fluids, and how far not. A closer 
examination shows us that for the interior of the fluid we can immediately apply our knowl- 
edge of the motion of a nonviscous fluid, but that care must be taken in considering the layers 
of the fluid in the immediate neighborhood of solid bodies. Frietion between fluid and solid 
body never comes into consideration in the fields of application to be treated here, because it 
is established by reliable experiments that fluids like water and air never slide on the surface 
of the body; what happens is, the final fluid layer immediately in contact with the body is 
attached to it (is at rest relative to it), and all the friction of fluids with solid bodies is therefore 
an internal friction of the fluid. Theory and experiment agree in indicating that the transition 
from the velocity of the body to that of the stream in such a case takes place in a thin layer of 
the fluid, which is so much the thinner, the less the viscosity. In this layer, which we call the 
boundary layer, the forces due to viscosity are of the same order of magnitude as the forces due 
to inertia, as may be seen without difficulty. It is therefore important to prove that, however 

` small the viscosity is, there are always in a boundary layer on the surface of the body forces 
due to viscosity (reckoned per unit volume) which are of the same order of magnitude as those 
due to inertia. Closer investigation concerning this shows that under certain conditions there 
may occur a reversal of flow in the boundary layer, and as a consequence a stopping of the fluid 
in the layer which is set in rotation by the viscous forces, so that, further on, the whole flow is 
changed owing to the formation of vortices. The analysis of the phenomena which lead to the 
formation of vortices shows that it takes place where the fluid experiences a retardation of flow 
along the body. The retardation in some cases must reach a certain finite amount so that a 
reverse flow arises. Such retardation of flow occurs regularly in the rear of blunt bodies; there- 
fore vortices are formed there very soon after the flow begins, and consequently the results 
which are furnished by the theory of nonviscous flow can not be applied. On the other hand, 
in the rear of very tapering bodies the retardations are often so small that there is no noticeable © 
formation of vortices. The principal successful results of hydrodynamics apply to this case. 
Since it is these tapering bodies which offer specially small resistance and which, therefore, 
have found special consideration in aeronautics under similar applications, the theory can be 
made useful exactly for those bodies which are of most technical interest. 


1 From this consideration one can calculate the approximate thickness of the boundary layer for each special case. 
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For the considerations which follow we obtain from what has gone before the result that 
in the interior of the fluid its viscosity, if it is small, has no essential influence, but that for 
layers of the fluid in immediate contact with solid bodies exceptions to the laws of a nonviscous 
fluid must be allowable. We shall try to formulate these exceptions so as to be, as far as possi- 
ble, in agreement with the facts of experiment. 


2. A further remark must be made concerning the effect of the compressibility of the 
fluid upon the character of the flow in the case of the motion of solid bodies in the fluid. All 
actual fluids are compressible. In order to compress a volume of air by 1 per cent, a pressure 
of about one one-hundredth of an atmosphere is needed. In the case of water, to produce an 
equal change in volume, a pressure of 200 atmospheres is required; the difference therefore is 
very great. With water it is nearly always allowable to neglect the changes in volume arising 
from the pressure differences due to the motions, and therefore to treat it as absolutely incom- 
pressible. But also in the case of motions in air we can ignore the compressibility so long as 
the pressure differences caused by the motion are sufficiently small. Consideration of compressi- 
bility in the mathematical treatment of flow phenomena introduces such great difficulties that 
we will quietly neglect volume changes of several per cent, and in the calculations air will be 
looked upon as incompressible. A compression of 3 per cent, for instance, oceurs in front of a 
body which is being moved with a velocity of about 80 m./sec. It is seen, then, that it appears 
allowable to neglect the compressibility in the ordinary applications to technical aeronautics. 
Only with the blades of the air screw do essentially greater velocities occur, and in this case the 
influence of the compressibility is to be expected and has already been observed. The motion 
of a body with great velocity has been investigated up to the present, only along general lines, : 
It appears that if the velocity of motion exceeds that of sound for the fluid, the phenomena are - 
changed entirely, but that up close to this velocity the flow is approximately of the same char- 
acter as in an incompressible fluid. 


3. We shall concern ourselves in what follows only with a nonviscous and incompressible 
fluid, about which we have learned that it will furnish an approximation sufficient for our 
applications, with the reservations made. Such & fluid is also called “the ideal fluid.” 

` What are the properties of such an ideal fluid? I do not consider it here my task to develop 
and to prove all of them, since the theorems of classical hydrodynamics are contained in all 
textbooks on the subject and may be studied there. I propose to state in what follows, for 
the benefit- of those readers who have not yet studied hydrodynamics, the most important 
principles and theorems which will be needed for further developments, in such a manner that 
these developments may be grasped. 1 ask these readers, therefore, simply to believe the 
theorems which I shall state until they have the time to study the subject in some textbook 
on hydrodynamics. 

The principal method of description of problems in hydrodynamics consists in expressing in 
formulas as functions of space and time the velocity of flow, given by its three rectangular com- 
ponents, u,v, w, and in addition the fluid pressure p. The condition of flow is evidently com- 
pletely known if u, v, w, and p are given as functions of x, y, z, and ¢, since then u, v, w, and p 
can be calculated for any arbitrarily selected point and for every instant of time. The direc- 
tion of flow is defined by the ratios of u, v, and w; the magnitude of the velocity is Vu? -v?+ wv. 
The “streamlines” will be obtained if lines are drawn which coincide with the direction of 
flow at all points where they touch, which can be accomplished mathematically by an inte- 
gration. If the flow described by the formulas is to be that caused by a definite body, then 
at those points in space, which at any instant form the surface of the body, the components of 
the fluid velocity normal to this surface must coincide with the corresponding components 
of the velocity of the body. In this way the condition is expressed that neither does the fluid 
penetrate into the body nor is there any gap between it-and the fluid. If the body is at rest 
in a stream, the normal components of the velocity at its surface must be zero; that is, the flow 
must be tangential to the surface, which in this case therefore is formed of stream lines. 
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4. In a stationary flow—that is, ina flow which does not change with the time, in which 
then every new fiuid particle, when it replaces another particle in front of it, assumes its velocity, 
both in magnitude and in direction and also the same pressure—there is, for the fluid particles 
lying on the same stream line, a very remarkable relation between the magnitude of the velocity, 
designated here by V, and the pressure, the so-called Bernouilli equation— 


p+ 5 V? = const. (1) 


(p is the density of the fluid, i. e., the mass of a unit volume). This relation is at once appli- 
cable to the case of a body-moving uniformly and in a straight line in a fluid at rest, for we are 
always at liberty to use for our discussions any reference system having a uniform motion in a 
straight line. If we make the velocity of the reference system coincide with that of the body, 
then the body is at rest with reference to it, and the flow around it is stationary. If now V 
is the velocity of the body relative to the stationary air, the latter will have in the new refer- 
ence system the velocity V upon the body (a man on an airplane in flight makes observations 
in terms of such a reference system, and feels the motion of flight as “wind”). 

The flow of incident air is divided at a blunt body, as shown in figure 1. At the point A 
the flow comes completely to rest, and then is again set in motion in opposite directions, tan- 
gential to the surface of the body: We learn from equation 
(1) that at such a point, which we shall call a “rest-point,”’ 


the pressure must be greater by 5 V? then in the undisturbed —s- 
fluid. We shall call the magnitude of this pressure, of which J) G 
we shall make frequent use, the “dynamical pressure,” and a-d f, 
shall designate it by g. An open end of a tube facing the 


stream produces a rest point of a similar kind, and there arises 
in the interior of the tube, as very careful experiments have a 
shown, the exact dynamical pressure, so that this principle ZI 


can be used for the measurement of the veloeity, and is in 
fact much used. The dynamical pressure is also well suited 
to express the laws of air resistance. It is known that this resistance is proportional to 


the square of the velocity and to the density of the medium; but g=5V ; so the law of air 
resistance may also be expressed by the formula 


W=c.F.g 2) 


where F is the area of the surface and e is a pure number. With this mode of expression it 
appears very clearly that the force called the “drag” is equal to surface times pressure differ- 
ence (the formula has the same form as the one for the piston force in a steam engine). This 
mode of stating the relation has been introduced in Germany and Austria and has proved use- 
ful. The air-resistance coefficients then become twice as large as the “absolute”. coefficients 
previously used. 

Since V? can not become less than zero, an increase of pressure greater than g can not, by 
equation (1), occur. For.diminution of pressure, however, no definite limit can be set. In 
the case of flow past convex surfaces marked increases of velocity of flow occur and in connection 
with them diminutions of pressure which frequently amount to 3q and more. 


Li 


Fic. 1.— Flow around a blunt body. 


5. A series of typical properties of motion of nonviscous fluids may be deduced in a useful 
manner from the following theorem, which is due to Lord Kelvin. Before the theorem itself 
is stated, two concepts must be defined. 1. The circulation: Consider the line integral of the 
‚velocity / V cos (V, ds). ds, which is formed exactly like the line integral of a force, which is 
called “the work of the force.” The amount of this line integral, taken over a path which 
returns on itself is called the circulation of the flow. 2. The fluid line: By this is meant a line 
which is always formed of the same fluid particles, which therefore shares-in the motion of the 
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fluid. The theorem of Lord Kelvin is: In a nonviscous fluid the circulation along every fluid 
line remains unchanged as time goes on. But the following must be added: 

(1) The case may arise that a fluid line is interseeted by a solid body moving in the fluid. 
If this occurs, the theorem ceases to apply. As an example I mention the case in which one 
pushes a flat plate into a fluid at rest, and then by means of the plate exerts a pressure on the 
fluid. By this a circulation arises which will remain if afterwards the plate is quickly withdrawn 
in its own plane. See figure 2. 

(2) In order that the theorem may apply, we must exclude mass forces of such a character 
that work is furnished by them along a path which returns on itself. Such forces do not ordi- 
narily arise and need not be taken into account here, where we are concerned regularly only 
with gravity. 

(3) The fluid must be homogeneous, i. e., of the same density at all points. We can easily 
see that in the case of nonuniform density circulation can arise of itself in the course of time ` 
if we think of the natural ascent of heated air in the midst of cold air.. The circulation increases 
continuously along a line which passes upward in the warm air and returns downward in the 
cold air. 

` Frequently the case arises that the fluid at the beginning is at rest or in absolutely uniform 
motion, so that the circulation for every imaginable closed line in the fluid is zero. Our theorem 
then says that for every closed line that can arise from one of the originally closed lines the 
circulation remains zero, in which we must make exception, as mentioned above, of those lines 
which are cut by bodies. If the line integral along every closed line is zero, the line integral 
for an open curve from a definite point O to an arbitrary point P is independent of the selection 
of the line along which the integral is taken (if this were not so, and if the 
integrals along two lines from O to P were different, it is evident that the 
1 line integral along the closed curve OPO would. not be zero, which contra- 

dicts our premise). The line integral along the line OP depends, therefore, - 
; _ since we will consider once for all the point 0 as a fixed one, only on the coordi- 
Pie ie br miodar nates of the point P, or, expressed differently, it is a function of these coor- 
tion and withdrawal of dinates. From analogy with corresponding considerations in the case of 
flat plate. fields of force, this line integral is called the ‘‘velocity potential,” and the 
particular kind of motion in which such a potential exists is called a “potential motion.” As 
follows immediately from the meaning of line integrals, the component of the velocity in a 
definite direction is the derivative of the potential in this direction. If the line-element is 
perpendicular to the resultant velocity, the increase of the potential equals zero, i. e., the sur- 
faces of constant potential are everywhere normal to the velocity of flow. The velocity itself 
is called the gradient of the potential. The velocity components u, v, w are connected with the 

potential by the following equations: 


ð$ 0b ð$ 
an? Say Y= de (3) 


The fact that the flow takes place without any change in volume is expressed by stating that 
as much flows out of every element of volume as flowsin. This leads to the equation 


az ay! az (4) 
In the case of potential flow we therefore have 

dp Fp $ 

tat rn (4a) 


as the condition for flow without change in volume. All functions ® (z, y, z, t), which satisfy 
this last equation, represent possible forms of flow. This representation of a flow is specially 
convenient for calculations, since by it the entire flow is given by means of the one function ®. 
The most valuable property of the representations is, though, that the sum of two, or of as 
many as one desires, functions ®, each of which satisfies equation (4a), also satisfies this equation, 
and therefore represents a possible type of flow (“superposition of flows”). 
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6. Another concept can be derived from the circulation, which is convenient for many con- 
siderations, viz, that of rotation. The component of the rotation with reference to any axis 
‘is obtained if the circulation is taken around an elementary surface of unit area in a plane 
perpendicular to the axis. Expressed more exactly, such a rotation component is the ratio of 
the circulation around the edge of any such infinitesimal surface to the area of the surface. The 
total rotation is a vector and is obtained from the rotation components for three mutually per- 
pendicular axes. In the case that the fluid rotates like a rigid body, the rotation thus defined 
comes out as twice the angular velocity of the rigid body. If we take a rectangular system of 
axes and consider the rotations with reference to the separate axes, we find that the rotation can 
also be expressed as the geometrical sum of the angular velocities with reference to the three axes. 

‘The statement that in the case of a potential motion the circulation is zero for every 
closed ‘fluid line can now be expressed by saying the rotation in it is always zero. The theorem 
that the circulation, if it is zero, remains zero under the conditions mentioned, can also now . 
be expressed by saying that, if these conditions are satisfied in a fluid in which there is no 
rotation, rotation can never arise. An irrotational fluid motion, therefore, always remains 
irrotational. In this, however, the following exceptions are to be noted: If the fluid is divided 
owing to bodies being present in it, the theorem under consideration does not apply to the 
fluid layer in which the divided flow reunites, not only in the case of figure 2 but also in the 
case of stationary phenomena as in figure 3, 
since in this case a closed fluid line drawn in 
front of the body can not be transformed into 
a fluid line that intersects the region where the 
fluid streams come together. Figure 3 shows 
‘four successive shapes of such a fluid line. This 
region is, besides, filled with fluid particles which 
have come very close to the body. We are 
therefore led to the conclusion from the stand- 
point of a fluid with very small but not entirely Ft ; 
vanishing viscosity that the appearance of vor- . FIG. 3.--Successive positions a line in flow around a solid 
tices at the points of reunion of the flow in the 
rear of the body does not contradict the laws of hydrodynamics. The three components of the 
rotation & n, ¢ are expressed as follows by means of the velocity components u, v, w 


ow òw _ ðu ðw , dv ðu 

£= Oy Oz” dz dz TS Sg ay 

If the velocity components are derived from a potential, as shown in equation (2), the rotation 
edb 8d 
azoy ðyðz 

7. Very remarkable theorems hold for the rotation, which were discovered by v. Helmholtz 
and stated in his famous work on vortex motions. Concerning the geometrical properties of the 
rotation the following must be said: 

At all points of the fluid where rotation exists the direction of the resultant rotation axes 
can be indicated, and lines can also be drawn whose directions coincide everywhere with these 
axes, just as the stream lines are drawn so as to coincide with the directions of the velocity. 
These lines will be called, following Helmholtz, “vortex lines.” The vortex lines through the 
points of a small closed curve form’a tube called a “vortex tube.” It is an immediate con- 
sequence of the geometrical idea of rotation as deduced above that through the entire extent 
of a vortex tube its strength—i. e., the circulation around the boundary of the tube—is constant. 
It is seen, in fact, that on geometrical grounds the space distribution of rotation quite inde- 
pendently of the special properties of the velocity field from which it is deduced is of the same 
nature as the space distribution of the velocities in an incompressible fluid. Consequently a 
vortex tube, just like a stream line in an incompressible fluid, can not end anywhere in the 
interior of the fluid; and the strength of the vortex, exactly like the quantity of fluid passing 
per second through the tube of stream lines, has at one and the same instant the same value 


(5; 


components, according to equation (5) vanish identically, since —— 
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throughout the vortex tube. If Lord Kelvin’s theorem is now applied to the closed fluid line 
which forms the edge of a small element of the surface of a vortex tube, the circulation along it 
is zero, since the surface inclosed is parallel to the rotation axis at that pomt. Since the circula- 
tion can not change with the time, it follows that the element of surface at all later times will 
also be part of the surface of a vortex tube. If we picture the entire bounding surface of a vortex 
tube as made up of such elementary surfaces, it is evident that, since as the motion continues 
this relation remains unchanged, the particles of the fluid which at any one time have formed 
the boundary of a vortex tube will continue to form its boundary. From the consideration 
of the circulation along a closed line inclosing the vortex tube, we see that this circulation—i.e., 
the strength of our vortex tube-—has the same value at all times. Thuis we have obtained the 
theorems of Helmholtz, which now can be expressed as follows, calling the contents of a vortex 
tube a “vortex filament”: “The particles of a fluid which at any instant belong to a vortex 
filament always remain in it; the strength of a vortex filament throughout its extent and for 
all time has the same value.” From this follows, among other things, that if a portion of the 
filament is stretched, say, to double its length, and thereby its cross section made one-half as 
great, then the rotation is doubled, because the strength of the vortex, the product of the rota- 
tion and the cross section, must remain the same. We arrive, therefore, at the result that the 
vector expressing the rotation is changed in magnitude and direction exactly as the distance 
between two neighboring particles on the axis of the filament is changed. 


8. From the way the strengths of vortices have been defined it follows for a space filled 
with any arbitrary vortex filaments, as a consequence of a known theorem of Stokes, that 
the circulation around any closed line is equal to the algebraic sum of the vortex strengths 
of all the filaments which cross a surface having the closed line as its boundary. If this closed 
line is in any way continuously changed so that filaments are thereby cut, then evidently the 
circulation is changed according to the extent of the strengths of the vortices which are cut. 
Conversely we may conclude from the circumstance that the circulation around a closed line 
(which naturally can not be a fluid line) is changed by a definite amount by a certain displace- 
ment, that by the displacement vortex strength of this amount will be cut, or expressed differ- 
ently, that the surface passed over by the closed line in its displacement is traversed by vortex 
filaments whose strengths add up algebraically to the amount of the change in the circulation. 

The theorems concerning vortex motion are specially important because in many cases 
it is easier to make a statement as to the shape of the vortex filaments than as to the shape of 
the stream lines, and because there is a mode of calculation by means of which the velocity 
at any point of the space may be determined from a knowledge of the distribution of the rota- 
tion. This formula, so important for us, must now be discussed. If T is the strength of a 
thin vortex filament and ds an element of its medial line, and if, further, 7 is the distance from 
the vortex element to a point P-at which the velocity is to be calculated, finally if a is the angle 
between ds and 7, then the amount of the velocity due to the vortex element is 

T ds sin a 

du=— Fa . (6) 
the direction of this contribution to the velocity is perpendicular to the plane of ds andr. The 
total velocity at the point P is obtained if the contributions of all the vortex elements present, 
in the space are added. The law for this calculation agrees then exactly with that of Biot- 
Savart, by the help of which the magnetic field due to an electric current is calculated. Vor- 
tex filaments correspond in it to the electric currents, and the vector of the velocity to the 
vector of the magnetic field. _ 

As an example we may take an infinitely long straight vortex filament. The contributions 
to the velocity at a point P are all in the same direction, and the total velocity can be deter- 
mined by a simple integration of equation (6). Therefore this velocity is 


+% . 
v r ds-sin a 
Ar r? 
© 
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As seen by figure 4, s=h ctg a, and by differentiation, ds = — de. Further r= na 
so. that 
da=— leos al) = (6a) 
=y; a sin ada = u zleos a Ich a) 


This result could be deduced in a simpler manner from the concept of circulation if we were 
to use the theorem, already proved, that the circulation for any closed line coincides with 
the vortex strength of the filaments which are inclosed by it. The circulation for every closed 
line which goes once around a single filament must therefore coincide with its strength. If 
the velocity at a point of a circle of radius A around our straight filament equals v then this 


circulation equals “path times velocity” =2rh.v, whence immediately follows v= ay The 


more exact investigation of this velocity field shows that for every point outside the filament 
(and the formula applies only to such points) the rotation is zero, so that in fact we are treat- 
ing the case of a velocity distribution in which only along the axis does rotation prevail, at 
all other points rotation is not present. 

For a finite portion of a straight vortex filament the preceding calculation gives the value 


v= gp(eos a, -- dos a) (6b) 


This formula may be applied only for a series of portions of vortices which together give an 
infinite or a closed line. The velocity field of a single portion of a filament would require 
rotation also outside the filament, in the sense that from the end of 
the portion of the filament vortex lines spread out in all the space 
and then all return together at the beginning of the portion. In the 
case of a line that has no ends this external rotation is removed, 
since one end always coincides with the beginning of another portion 
of equal strength, and rotation is present only where it is predicated Fis. 4.—Velocity-field due to infinite 
in the calculation. rectilinear vortex, 


9, If one wishes to represent the flow around solid bodies in a fluid, one can in many cases 
proceed by imagining the place of the solid bodies taken by the fluid, in the interior of which 
disturbances of flow (singularities) are introduced, by which the flow is so altered that the 
boundaries of the bodies become streamline surfaces. For such hypothetical constructions 
in the interior of the space actually occupied by the body, one can assume, for instance, any 
suitably selected vortices, which, however, since they are only imaginary, need not obey the 
laws of Helmholtz. As we shall see later, such imaginary vortices can be the seat of lifting 
forces. Sources and sinks also, i. e., points where fluid continuously appears, or disappears, 
offer a useful method for cohstructions of this kind. While vortex filaments can actually 
occur in the fluid, such sources and sinks may be assumed only in that part of the space which 
actually is occupied by the body, since they represent a phenomenon which can not be realized. 
A contradiction of the law of the conservation of matter is avoided, however, if there are assumed 
to be inside the body both sources and sinks, of equal strengths, so that the fluid produced by 
the sources is taken back again by the sinks. 

The method of sources and sinks will be described in greater detail when certain practical 
problems are discussed; but at this point, to make the matter clearer, the distribution of veloci- 
ties in the case of a source may be described. It is very simple, the flow takes place out from 
the source uniformly on. all sides in the direction of the radii. Let us describe around the point 
source a concentric spherical surface, then, if the fluid output per second is Q, the velocity at 
the surface is 


P= yy 0 
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the velocity therefore decreases inversely proportional to the square of the distance. The 
flow is a potential one, the potential comes out (as line-integral along the radius) 


Q 


$= const. — Arr (7a) 


If a uniform velocity toward the right of the whole fluid mass is superimposed on this 
velocity distribution—while the point source remains stationary—then a flow is obtained 
which, at a considerable distance from the source, is in straight lines from left to right. The 
fluid coming out of the source is therefore pressed toward the right (see fig. 5); it fills, at some 
distance from the source, a cylinder whose diameter may be determined easily. If V is the 
velocity of the uniform flow, the radius r of the cylinder is given by the condition Q= ar’. V. 
All that is necessary now is to assume on the axis of the source further to the right a sink of 
the same strength as the source for the whole mass of fluid from the source to vanish in this, 
and the flow closes up behind the sink again exactly as it opened out in front of the source. 
In this way we obtain the flow around an elongated body with blunt ends. 


10. The special case when in a fluid flow the phenomena in all planes which are parallel 
to a given plane coincide absolutely plays an important rôle both practically and theoretically. 
If the lines which connect the corresponding 
points of the different planes are perpendicular to 
the planes, and all the streamlines are plane 
curves which lie entirely in one of those planes, 
we speak of a uniplanar flow. ‘The flow around 
a strut whose axis is perpendicular to the direc- 
tion of the wind is an example of such a motion. 

The mathematical treatment of plane poten- 
tial flow of the ideal fluid has been worked out 
specially completely more than any other prob- 
lem in hydrodynamics. This is due to the fact 
that with the help of the complex quantities 
(e+iy, where i=-/—1, is called ‚the imaginary 

5 iti A axe. Unit) there can be deduced from every analytic 
Fig. 5.—Superposition of uniform flow and that caused by a soure: function a case of flow of this type which is incom- 
pressible and irrotational. Every real function, ® (x, y) and'Y (x, y), which satisfies the relation 


b+iV=fletiy), (8) 
where f is any analytic function, is the potential of such a flow. This can be seen from these 
considerations: Let «-+iy be put=z, where z is now a “complex number.” Differentiate equa- 
tion (8) first with reference to z and then with reference to y, thus giving 


ô 3x dfdz_df 
da da dzdx de 
db, av _dfaz _ df _,0®_av 
dy ay dzdy de ‘dx dy 
In these the real parts on the two sides of the equations must be equal and the imaginary 
parts also. If @ is selected as the potential, the velocity components u and v are given by 
_9b_ ay db _ a 
“= dy? dy de 
ðu 
ðc 


(9) 


If now we write the expressions 


+2 (continuity) and 2g (rotation) first in terms of 
® and then of F, they become 


APPLICATIONS OF MODERN HYDRODYNAMICS TO AERONAUTICS. 


du , 3v _ FS oe or ay 

Oz Oy ðL dy? dydx dxdy 
=0 

av ðu _ Pb Ph OY HY 

Ox Oy dydx dxdYy ax? 8y? 


(10) 


=0 
It is seen therefore that not only is the motion irrotational (as is self-evident since there is 
a potential), but it is also continuous. The relation +5 7° besides corresponds exactly 


to our equation (4a). Since it is satisfied also by Y, this can also be used as potential. 
The function Y, however, has, with reference to the flow deduced by using © as potential, a 
special individual meaning. From equation (8) we can easily deduce that the lines Y= const. 
are parallel to the velocity; therefore, in other words, they are streamlines. In fact if we put 


av 

dy _ dz v 

adv Far +2 mel 0, then a ove m 
oy 


- which expresses the fact of parallelism. The lines Y= const. are therefore perpendicular to the 
lines ®=const. If we draw families of lines, ®=const. and Y=const. for values of and Y 
which differ from each other by the same small amount, it follows from the easily derived 


equation db +idv = df 7, (de + idy) that the two bundles form a square network; from which fol-: 


lows that the a curves of the network again form an orthogonal and in fact a square 
network. This fact can be used practically in drawing such families of curves, because an error 
in the drawing can be recognized by the eye in the wrong shape of the network. of diagonal 
curves and so can be improved. With a little practice fairly good accuracy may be obtained 
by simply using the eye. Naturally there are also mathematical methods for further improve- 
ment of such networks of curves. The function Y, which is called the “stream function,” 
has another special meaning. If we consider two streemlines ¥=W, and ¥=W,, the quantity 
of fluid which flows between the two streamlines in a unit of time in a region of uniplanar flow 
of thickness 1 equals ¥, — Y. In fact if we consider the flow through a plane perpendicular to 


the X-axis, this quantity is 
o- ar (a [an 


The numerical value of the sfream function coincides therefore with the quantity of fluid which 
flows between the point x, y and the streamline V=o0. 
As an example let the function , 
$ +it = A(z +iy)” 
be discussed briefly. It is simplest in general to ask first about the streamline ¥=o. As is 
well known, if a transformation is made from rectangular coordinates to polar ones r, y, (x-+ty)* 
=r? (cos np+i sin ne). The imaginary part of this expression is ir” sin ne. This is to be 
put equal tow. W=o therefore gives sin ng =o, i. e., np=0, n,2r,etc. The streamlines ¥=0 


are therefore straight lines through the origin of coordinates, which make an angle ane with 


each other, the flow is therefore the potential flow between two plane walls making the angle 
a with each other. The other streamlines satisfy the equation r" sin ne=const. The veloci- 
ties can be obtained by differentiation, e. g., with reference to x: 


oğ 


+ =u-Ww=An (x +) na = An (cos (n-1) y +å sin (n-1)o} 


0% 
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For r=o this expression becomes zero or infinite, according as n is greater or less than 11, i. e., 
according as the angle a is less or greater than r(=180°). Figures 6 and 7 give the streamlines 


for a = =45° and ar =270°, corresponding to n=4 and S. In the case of figure 7 the velocity, 


as just explained, becomes infinite at the corner. It would be expected that in the case of 
the actual flow some effect due to friction would enter. In fact there are observed at such 
corners, at the beginning of the motion, great velocities, and immediately thereafter the for- 
mation of vortices, by which the motion is so changed that the velocity 
at the corner becomes finite. 


It must also be noted that with an equatiorf 
( p+iq=o(e+ry) (11) 
the z-y plane can be mapped upon the p-g plane, since to every pair 
iN 4 


of values z,y a pair of values p,q corresponds, to every point of the x-y 
plane corresponds a point of the p-g plane, and therefore also to every 
Fic. t.-Unipianar flow be- element of a line or to every curve in the former plane a linear element 
tween plane walls making : i fi . 
an angle a=45° with each 22d & curve in the latter plane. The transformation keeps all angles 
other. unchanged, i. e., corresponding lines intersect in both figures at the same 
angle. 


By inverting the function p of equation (11) we can write 


t+y=x (p +i) 
and therefore deduce from equation (8) that 


$ +it =f {x(pt+ig]=F (p+) (12) 
® and Y are connected therefore with p and g by an equation of the type of equation (8), and 
hence, in the p-q plane, are potential and stream functions of a flow, and further of that, flow 
which arises from the transformation of the 6, Y network in the z-y plang into the p-g plane. 

This is a powerful method used to obtain by transformation from a known simple flow 
new types of flow for other given boundaries. Applications of this will be given in section 14. 


11. The discussion of the principles of the hydrodynamics of nonviscous fluids to be 
applied by us may be stopped here. I add but one considera- 
tion, which has reference to a very useful theorem for obtaining 


the forces in fluid motion, namely the so-called “momentum theo- de 
rem for stationary motions.” 
We have to apply to fluid motion the theorem of general 
mechanics, which states that the rate of change with the time 
of the linear momentum is equal to the resultant of all the ex- 
ternal forces. To do this, consider a definite portion of the 
i 


fluid separated from the rest of the fluid by a closed surface. 

This surface may, in accordance with the spirit of the theorem, 

be considered as a “fluid surface,” i. e., made up always of Fie. 7—Uniplanar flow around plane 
the same fluid particles. We must now state in a formula the walls making an angle 270° with each 
change of the momentum of the fluid within the surface. If, as ics 

we shall assume, the flow is stationary, then after a time dt every fluid particle in the interior 
wili be replaced by another, which has the same velocity as had the former. On the boundary, 
however, owing to its displacement, mass will pass out at the side where the fluid is approaching, 
and a corresponding mass will enter on the side away from which the flow takes place. If dS 
is the area of an element of surface, and v, the component of the velocity in the direction of 
the outward drawn normal at this element, then at this point dm=pdS . v, dt. If we wish 
to derive the component of the “impulse”’—defined as the time rate of the change of momen- 
tum—for any direction s, the contribution to it of the element of surface is 


d 
dJ, =v = PS < Unds (13) 
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With this formula we have made the transition from the fluid surface to a corresponding solid 
“control surface.” 

The external forces are compounded of the fluid pressures on the control surface and the 
forces which are exercised on the fluid by any solid bodies which may be inside of the control 
surface. If we call the latter P, we obtain the equation 


zP.=ffp . cos (n, $) . dS +p f fvnvdsS (14) 


for the s component of the momentum theorem. The surface integrals are to be taken over the 
entire closed control surface. The impulse integral can be limited to the exit side, if for every 
velocity v, on that side the velocity v,’ is known with which the same particle arrives at the 
approach side. Then in equation (13) dJ is to bé replaced by 


ER, on de (138) 


The applications given in Part. II will furnish illustrations of the theorem. 
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"APPLICATIONS OF MODERN HYDRODYNAMICS TO AERONAUTICS. 


By L. PRANDTL. 


PART II. 
APPLICATIONS. 


A. DISTRIBUTION OF PRESSURE ON AIRSHIP BODIES. 


12. The first application of hydrodynamical theory to be tested by experiment in the 
Göttingen Laboratory referred to the distribution of pressure over the surface of models of 
airships. We can construct mathematically the flow for any number of varieties of sectional 
forms of bodies of revo'ution of this kind if we place along an axis parallel with the direction 
of the air current any suitable distribution of sources and sinks, taking care that the total 
strength of the sources and sinks are the same. According to the intensity of the uniform 
motion which is superimposed upon the flow from the sources, we obtain from the same system 
of sources and en bodies of different thicknesses. In order to obtain the smoothest possible 
shapes, the sources and sinks are generally distributed contin- 
uously along the axis, although single-point sources are allowable. 

` In the case of continuously distributed’sources and sinks the 
method of procedure is briefly this: The abscissas of the single 
sources are denoted by &, the intensity of the source per unit of 
length by /(&), in which positive values of f(£) denote sources, nega- 
tive values sinks. The condition that makes the stream from the — 
sources self-contained is expressed by the equation 


Fig. 8,--Explanation of quantities used 
in calculation of streamline shapes. 


l 
Abseissae: Position of source. J Sf®@dt=o. 
o 


Ordinates: Intensity of source per unit 
length, i. e.f (£). s R f š 
By simply adding the potentials due to the single elementary 


sources f(E) dé, i. e., in this case by integrating them, the total flow due to the sources will be 
given by the potential defined by the following formula 


i 
1 d (15 


in which r= -+/(@—#)?+ 4’, and y is the perpendicular distance from the axis of the point for 
which the potential is calculated, x is the abscissa along the axis measured from the same 
origin as £. (See fig. 8.) There must be added to this potential that due to the uniform flow 
with the velocity V, viz, ®,= Vz. The total potential is then =, +,; and therefore the veloc- 
ity parallel to the axis is u= a? = V+ 5 and the sidewise (radial) velocity is v= 2, 7 
In order to calculate the streamlines one could perform an integration of the direction 
given by u and v. These lines are obtained more conveniently, in this case also, by means of 
the stream function. (See sec. 10.) In the case of flow symmetrical with reference to the axis, 
such as is here discussed, one can take as stream function the quantity of fluid flowing inside 
the circle drawn through the point z, y in a plane perpendicular to the axis and having its center 
172 
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on the axis. The amount of fluid delivered by the sources which lie up the stream is purposely 
deducted from this. It is not difficult to see that all points of the X-Y plane, through whose 
parallel circles the same amount of fluid flows per second—after deduction of the sources— 


must lie on one and the same streamline, for evidently there is no flow, either in or out, through | 


the surface formed by the streamlines drawn through the points of any one parallel circle 
(since the flow is along the surface); therefore the quantity of fluid flowing within this surface is 
constant, so far as it is not increased by the sources. From the meaning of the stream function, 
to determine which the velocity must be integrated over a surface, it follows that the stream 
function of a flow due to two or more causes is at every point the sum of the stream functions 
of the several partial flows. For a continuous distribution of sources therefore the stream 
function y is. obtained by an integration exactly as was the potential. According to our 
premise the surface of the body is designated simply by the value y=o. The formulas are 
obtained as follows: 

The flow from a simple source through a circle passing through a point lying to the right of 
the source is, writing r= V2? +4, 


,_ [p _ [?Qx _ Qu (rydy_Q/ 2) 
Y = [uorudy= fase 2rydy= 5 Í n = (: z 
From this, in accordance with what has been said, the quantity Q must be subtracted, so that 
Q/ 2) ; 
d = 
y=w'—9=—S(145 (16) 
For points lying to the left of the source we obtain from the integral i 


(44) 


which coincides with formula (16); this holds, then, everywhere. 
For the assumed continuous distribution of sources we obtain 


1 = 
n=- [jo (14778) ae (17) 

o r 
in which r= -/(c—?+y?. To this stream function of the sources must now be added that due 
to the paral!l flow Y, = Vm (18) 


Putting the total stream function Y, +Y, =Y equal to zero, gives the equation of the surface 

“of the body around which the flow takes place. Putting ¥;+¥,= C gives any other streamline. 
It is evident that, with the same distributian of sources, a wholè group 
of body surfaces can be obtained, depending upon the choice of the ratio 
of the intensity of the sources to the strength of the parallel flow. 

The determination is best made practically by graphical methods, for 
‘instance, by laying off the curves x= const. in a system of coordinates 
consisting of y and —W, which can be obtained at once from a calculation 
by tables for the stream function ¥,. If we intersect these curves by 
parabolas corresponding to the equation —Y= Vry — C, we obtain at 

_ once a contour (for C=0), or some external or internal streamline (for 
C>o or C<o). The parabola may be drawn upon transparent paper, and Fic. 9.—The curves are for 
then by displacing the parabola along the ¥ axis we can at once obtain different values of z= 
from figure 9 the values of y corresponding to any x. pi 

In this manner a former colleague of mine,-who unfortunately fell immediately at the 
beginning of the war, Dr. G. Fuhrmann, calculated the shapes of bodies corresponding to a 
series of source distributions, and on the one hand he determined the distribution of pressure 
over the surface of these bodies by means of the Bernouilli equation (see sec. 4) l 


p=pot {V> Ut) ) (19) 


ı The velocities u and v may be obtained from the potential, but also from the stream function y; för u= z = and v= l 8%, 
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and on the other he constructed models according to these drawings and measured the pressure 
distribution over them when placed in a wind tunnel. The agreement was altogether surpris- 
ingly good, and this success gave us the stimulus to seek further relations between theoretical 
hydronamics and practical aeronautics. The work of Fuhrmann was published in Jahr-b. 
der Motorluftschiff-Studien Gesellsch., Volume V, 1911-12 (Springer, Berlin), and contains 
a large number of illustrations. Four of the models investigated are shown here. The upper 
halves of figures 10 to 13 show the streamlines for a reference system at rest with reference to 
the undisturbed air, the lower halves the streamlines for a reference system attached to the 
body.. The distribution of the source intensities is indicated on the axis. The pressure dis- 
tributions are shown in figures 14 to 17. The calculated pressure disfributions are indicated 
by the lines which are drawn full, the individual observed pressures by tiny circles? _ . 

It is seen that the agreement is very complete; at the rear end, however, there appears a 
characteristic deviation in all cases, since the theoretical pressure distribution reaches the full 
dynamical pressure at the point where the flow reunites again, while actually this rise in pres- 
sure, owing to the influence of the layer of air retarded by friction, remains close to the surface. 

As is well known there is no resistance for the theoretical flow in a nonviscous fluid. The 
actual drag consists of two parts, one resulting from all the normal forces (pressures) acting 
on the surface of the body, the other from all the tangential forces (frietion). The pressure 
resistance, which in this case can be obtained by integration of the pressure distribution over 
the surface of the body, arises in the main from the deviation mentioned at the rear end, and 
is, as is known, very small. Fuhrmann’s calculations gave for these resistances a coefficient, 
with reference to the volume of the body, as shown in the following table: € 


| Models esse nen en ee I Ir IM Iv | 


i ois NAINE ISEE E S N T ETEN OEA TA = 0. 0170 0. 0123 0. 0131 0.0145 


This coefficient is obtained from the following formula: 
Drag w=k, U7" q 
where U designates the volume and g the dynamical pressure. 
The total resistance (drag) was obtained for the four models by means of the balance; 


the difference between the two quantities then furnishes the frictional resistance. The total 
drag coefficients were: ? 


Model cerre oyin see see Seevs essen. I I III Iv 


ids ne a desta es Nol ests volar eae =| 0.0840 | 0.0220 | 0.0246 | 0.0248 


With greater values of VZ than were then available for us, the resistance coefficients 
become nearly 30 per cent smaller. For purposes of comparison with other cases it may be 
mentioned that the “maximum section” was about 2/5 of U?®, The surface was about seven 
times U?®; from which can be deduced that the total resistance of the good models was not 


- greater than the friction of a plane surface having the same area. The theoretical theorem 


that in the ideal fluid the resistance is zero receives in this a brilliant confirmation by experiment. 
B. THEORY OF LIFT. 


13. The phenomena which give rise to the lift of an aerofoil may be studied in the simplest 
manner in the case of uniplanar motion. (See sec. 10.) Such a uniplanar flow would be ex- 
pected obviously in the case that the wing was unlimited at the sides, therefore was ‘‘infinitely 


2 In the wind tunnel there was a small pressure drop in the direction ofits length. In order to eliminate the effect of this, the pressures toward 
the fore had to be diminished somewhat and those aft somewhat increased. - 
` 3 After deduction of the horizontal buoyancy. 
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Fig. 10. Fig. 11. 
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nn ay a es 
FIG. 12. Fig. 13, 


Four airship models as derived by Fuhrmann by combination of sources and uniform fiow. Distribution of sources indicated on axis. Upper half: 
Streamlines relative to undisturbed air. Lower half: Streamlines relative to airship. 


Fia. 16. Fig. 17. 


Pressure distribution over airships of figures 10 to 13. Full lines represent calculated values; small circles, points as found by observation in a 
wind-tunnel. 
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long,” and throughout exhibited the same profile and the same angle of attack. In this case all 
the sections will be alike in all respects and each one can be considered as a plane of symmetry. 
The infinitely long wing plays an important part therefore in the considerations of the theoretical 
student. It is not possible to realize it in free air, and marked deviations from the infinitely 


‘long wing are shown even with very long wings, e. g., those having an aspect ratio of 1:10. In 


laboratories, however, the infinitely long wing, or uniplanar flow, may be secured with good 
approximation, if a wing having a constant profile is placed between plane walls in a wind 
tunnel, the walls running the full height of the air stream. In this case the wing must extend 
close to the walls; there must be no gap through which a sensible amount of air can flow. We will 
now discuss such experiments, and first we shall state the funda- 
mental theory of uniplanar flow. 

Since, as explained in section 4, in a previously undisturbed 
fluid flow, the sum of the static and dynamic pressures is con- 


stant: p+5V?=const., in order to produce lift, for which the 


pressure below the surface must be increased and that above 
diminished, such arrangements must be made as will diminish 
the velocity below the wing and increase it above. The other 
method of. producing such pressure differences, namely, by 
Fig. 18.—Deduction of the Kutta formula, causing a vortex region above the surface placed hie a kite 
‘ihiplanar flow aronind ninlte wing oblique to the wind, by which a suction is produced, does not 
come under discussion in practical aeronautics owing to the great resistanceit setsup. . Lanchester 
has already called attention to the fact that this lifting current around the wing arises if there is 
superimposed upon a simple potential flow a circulating flow which on the pressure side runs 
against the main current and on the suction side. with it. Kutta (1902) and Joukowski (1906) 
proved, independently of each other, the theorem that the lift for the length Z of the wing is 


=r V1 (20) 


in which T is the circulation of the P flow. It may be concluded from this formula 
that in a steady fluid flow lift is not possible unless there is motion giving rise to a circulation. 
In uniplanar flow in an ideal fluid this lift does not entail 
any drag. 

The proof of the Kutta-Joukowski formula is generally NN 
deduced by applying the momentum theorem to a circular Sf 
cylinder of large radius whose.axis is the medial line of the 
wing. The circulatory motion, which could be obtained 


numerically close to the wing only by elaborate mathematical — L—— 


processes, is reduced at a great distance from the wing to a a... 
motion which agrees exactly with the flow around a rectilinear a a eee a 


vortex filament (see sec. 8), in which, therefore, the single Fre. 1».—Uniplanar uniform flow around eir« 
particles describe concentric ons The velocity around a NE 


circle of radius 2 is, then, Per For an element of surface J. Rd@ (see fig. 18) the normal 


component of the velocity is V cos 6, the mass flowing through per second dm=plRV cos bdo. 


- If we wish to apply the momentum theorem for the vertical components, i. e., those perpen- 


dicular to the direction of V, then this component of the velocity through the element of surface 
must be taken. This, obviously, is v cos 6, taken positive if directed downward; the total 
impulse, then, is 


J= fe cos 0dm = pl R Vv f 5 cos? 6dé. 


The integral equals r, and therefore introducing the value of v 


J=5 oVTl. 
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Since the resulting impulse is directed downward (the upward velocity in front of the wing is 
changed into a downward one behind the wing), this means that the reaction of the fiuid against 
the wing is a lift of the wing upward. The amount of the impulse furnishes, as is seen by re- 
ferring to formula (20), only half the lift. The other half comes from the pressure differences 
on the control surfaces. Since, for a sufficiently large R, v can always be considered small com- 


pared with V, neglecting 5 v, the pressure p is given, according to the ernouilli equation, by 
p=pot 8-8 (V40 sin +9? cos? }=p.— nV sin 6. 


A component of this, obtained by multiplying by sin 6, acts vertically on the surface element 
IRdo. The resulting force Dis, then, 3 

D=plRVo |T sin? ads. 
This integral also equals 7, so that here also 


= 5 PVT 
its direction is vertically up. The total lift, then, is 
A=J+D =p Vr. 


14. For the more accurate analysis of the flow around wings the complex functions (see 
sec. 10) have been applied with great success, following the procedure of Kutta. Very different 


Q 


Fic. 20.—Uniplanar flow around circular Fig. 21.—Superposition of two preceding 
cylinder considered as a columnar vortex flows. 
of strength F. 


methods have been used. Here we shall calculate only one specially simple case, in which the 
flow will be deduced first around a circular cylinder and then calculated for a wing profile by a 
transformation of the circular cylinder and its flow, using complex functions. 

The flow around a circular cylinder has long been known. If the coordinates in the plane 
of the circle are p and q, and if we write p+iq=t, the potential and stream functions for the 
ordinary symmetrical flow around the circular cylinder are given by the very simple formula 


2,4i0,=V(t+5) 5 (21) 


It is easily seen by passing to polar coordinates that, for r=a, ¥,=0, and that therefore the 
circle of radius a is a streamline. Further, for the p axis, y, =0, i. e., this is also a streamline. 
The whole flow is that shown in figure 19. To this flow must be added the circulation flow 
expressed by the formula 


. _ 49. 
t= 5° log ¢ ° (22) 
which, as shown in figure 20, is simply a flow in concentric circles with the velocity a The 


combination of the two flows, i. e., the flow for the sum of the expressions in equations (21) 
and (22), is shown in figure 21. It is seen that the rest point is moved down an amount de. 
By.a suitable choice of the circulation this can be brought to any desired point. 


sat(e) cos d-+i (-#) sin 3 6 į log t=—8+i log r. 
20167—23-——12 
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We must now discuss the transformation of this flow to a wing profile. For this purpose 
manifold means are possible. The simplest is furnished by a transformation according to the 


equation 


eantiy-1+ 


2 
a8 


By this the circle of diameter AB=2b in the ¢ plane (as we shall for brevity’s sake call the p, q 
plane) is transformed into a straight line A’ B’ of the length 4b along the X axis, and concen- 


180°; 


Fic. 22. 


2 
Conformal transformation of z plane into ¢ plane by ante 


FIG, 23. 


tric circles around the former become ellipses, the radii become hyperbolas. All the ellipses 
and hyperbolas have their foci at the ends of the straight line, this forming a confocal system. 
Figures 22 and 23 illustrate the transformation. It may be mentioned, in addition, that the 
interior of the circle in figure.22 corresponds to a continuation of the meshwork in figure 23 
through the slit A’ B’, whose form agrees with the meshwork as drawn. Any circle through 


Fie, 24.—Illustrations of Joukowski sections. 


the points AB is thereby transformed into an arc of 
a circle passed over twice, having an angle subtended 


_ at the center equal to 4. 


Many different results may now be obtained by 
means of this mapping, according to the position 
which the circle, around which the flow takes place 
according to equations (21) and (22), bears to the di- 
ameter AB of the circle of figure 22. If the diameter 
AB is made to coincide with any oblique diameter of 
the circular section of the cylinder, we obtain a flow 
around an oblique plate whose angle of attack coin- 
cides with the inclination of the line AB. If the di- 


ameter AB is selected somewhat smaller, so that both points lie inside the circle symmetrically 
on the diameter, the flow around ellipses is obtained. If, however, the diameter AB coincides 
with a chord of the circle around which the original flow was, which, for example, may lie below 
the center, the flow around a curved plate forming an arc of a circle is obtained. By selection 
of various points in the interior of the original circle forms of diverse shapes are obtained. The 
recognition of the fact that among these forms very beautiful winglike profiles may be found we 
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owe to Joukowski. These are obtained if the point B is selected on the boundary of the original 
circle and the point A inside, and somewhat below the diameter through the point B. Figure 
24 gives illustrations of such Joukowski profiles. 

In order that the flow may be like the actual one, in the cases mentioned the circulation 
must always be so chosen that the rear rest 


point coincides with the ‚point B, or, r- — > —_—— Z~ 


spectively, with the point on the original 


circle which lies nearest this point. In en 
this case there will be, after mapping on EEE oa O 


the z plane, a smooth flow away from the cl 2~~ 
trailing edge, as is observed in practice. It 
is therefore seen that the circulation must ean he Wiig goa as, ee 


be taken greater according as the angle of i a ee TEERSETE 
attack is greater, which agrees with the ob- DO  SEESGHEEI 


servation that the lift increases with increas- SS ee EEE 
ing angle of attack. 


The transformation of the flows shown 
in figures 19 to 21 into wing profiles gives 
illustrations of streamtines as shown in figures 25 to 27—figure 25, simple potential flow; 
figure 26, circulation flow; figure 27, the actual flow around a wing obtained by superposition 
of the two previous flows. 

We are, accordingly, by the help of such constructions, in the position of being able to 

calculate the velocity at every point in the 
neighborhood of the wing profile, and with it 


the pressure. In particular, the distribution 
of pressure over the wing itself may be cal- 
culated. 
MUU ay My assistant, Dr. A. Betz, in the year 
1914 worked out the pressure distribution for 
l a Joukowski wing profile, for a series of angles 
of attack, and then in a wind tunnel meas- 
ured the pressure distribution on a hollow 
model of such a wing made of sheet metal, 


side walls of the height of the tunnel being 

Fic. 26..- Transformation of circulatory flow, figure 20. introduced so as to secure uniplanar flow 

The results of the measurements agreed in a very satisfactory manner with the calculations, 

only—as could be well explained as due to friction—the actual circulation was always slightly 

less than that calculated for the same angle of attack. If the pressure distributions would 
be compared, not for the same angles of 


attack, but for the same amouht of circu- IT — In I 
lation, the agreement would be noticeably SS — 
better. The pressure distributions ‘are ——— III 
‘shown in figures 28 to 30, in which again > 
the full curves correspond to the measure- TE Woods, 


ments and the dashes to the calculated Lye 
pressures. Lift and drag for the, wing = 

were also obtained by the wind-tunnel bal- ann 
ance. In order to do this, the middle part ee PIERRE 
of the wing was isolated from the side parts, Zoe 
which were fastened to the walls of the tun- 

nel by carefully designed labyrinths, so that 
within a small range it could move without 
friction. The result of the experiment is shown in figure 31.. The theoretical drag is zero, that 


obtained by measurement is very small for that region where the wing is “good,” but sensibly 
larger for too large and too small angles of attack. The lift is correspondingly in agreement 


Fig. 25.— Transformation of simple potential flow, figure 19. 


Frc. 27.—Transformation of superposition of the two flows, figure 21. 
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with the theoretical value in the good region, only everywhere somewhat less. The deviations _ 
of drag as well as of lift are to be explained by the influence of the viscosity of the fluid. The 
agreement on the whole is as good as can be expected from a theory which neglects completely 
the viscosity. 

For the connection between the angle of incidence «œ and the circulation which results from 
the condition discussed above calculations give the following result for the lift: 

(1) The Kutta theory gives for the thin plane 
plate the formula 


A=bt.rp V sin a (23) 
The lift coefficient C, is defined by the equation 


Oe where 1-5 pV? 


and therefore 
(,=2r sin a (24) 


(2) For the circularly curved plates having an 
angle of arc 48 subtended at the center (see figure 23) 
we have, according to Kutta, if a is the angle of 
attack of the chord, 


(25) 


which, for small curvatures, becomes 27 sin (a+8); 
this can be expressed by saying that the lift of the 
circularly curved plate is the same as that of a plane 
REN which touches the former at a point three-fourths of 
the distance around the arc from its leading edge. 
For the Joukowski profiles and for others the 
‘formulas are less simple. v. Mises showed in 1917 
that the increase of C, with the angle of attack, i. e., 


‘a 


2 » is greater for all other profiles than for the flat 


plate, and is the greater the thicker the profile. But 
the differences are not marked for the profiles occur- 
ring in practice. 

The movement of the center of pressure has 
also been investigated theoretically. With the 
plane plate, in the region of small angles, it always 
lies at one-fourth of the width of the plate; with 
circularly curved thin plates its position for small 
angles is given by the following law: 


t tan æ 
F16. 30. %= 7 tan a+ tan ĝ (26) 


Pressure distribution over a Joukowski wing, different anglesof . m = . i 
attack. Full lines give results of wind-tunnel tests; dashed In which żis the chord of the plate, and x, is the dis- 


lines, calculated values. tance measured from the center of the plate. The 
fact that the movement of the center of pressure in the case of “good” angles of attack of the 
profiles agrees with theory is proved by the agreement of the actual pressure distribution with 
that calculated. In the case of thin plates a less satisfactory agreement as respects pressure 
distribution is to be expected because with them in practice there is a formation of vortices at 
the sharp leading edges, while theory must assume a smooth flow at this edge. 


15. That a circulatory motion is essential for the production of lift of an aerofoil is defi- 
nitely established. The question then is how to reconcile this fact with the proposition that 
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the circulation around a fluid line in a nonviscous fluid remains constant. If, before the motion 
begins, we draw a closed line around the wing, then, so long as everything is at rest, the circula- 
tion certainly is zero. Even when the motion begins, it can not change for this line. The ex- 
planation of why, in spite of this, the wing gains circulation is this: At the first moment of the 
motion there is still no circulation present, the motion takes place approximately according to 
figure 25, there is a flow at high velocity around the 
trailing edge. (Seesec.10.) This motion can not, how- 
ever, continue; there is instantly formed at the trailing 
edge a vortex of increasing intensity, which, in accord- 
ance with the Helmholtz theorem that the vortex is 
always made up of the same fluid particles, remains with 
the fluid asit passes on. (See fig. 32.) The circulation 
around the wing and vortex, taken together, remains 
equal to zero; there remains then around the wing a cir- 
culation equal and opposite to that of the vortex which 
has gone off with the current. Therefore vortices will 
be given off until the circulation around the wing is of ` 
such a strength as to make the fluid fiow off smoothly 
from the trailing edge. If by some alteration of the 
angle of attack the condition for smooth flow is ‘dis- 
turbed, vortices are again given off until the circulation 
reaches its new value.. These phenomena are com- 
pleted in a comparatively short distance, so the full lift 
is developed very quickly. 

In the pictures of flow around a wing, e. g., figure 27, 
one sees that the air in front of the wing flows upward 
against the reaction of the lift. The consideration of 
momentum has shown that half of the impulse is due to 
the oncoming ascending current. This fact needs some 
further explanation. The best answer is that given by 
Lanchester,’ who shows that for the production of lift Fie. 31.—Values of lift and drag coeflieientsofa Toukowski 
the air mass at any time below the wing must be given wing as obtained in wind-tunnel tests and hy theory. 
an acceleration downward. The question he asks is: What kind of a motion arises if for a short 
time the air below the wing is accelerated downward, then the wing is moved forward a bit 
without pressure, then the air is again accelerated, and so on? The space distribution of the 
accelerations is known for the case of a plane plate, infinitely extended at the sides, accelerated 
l P et from rest;. the pattern of the accelera- 
7 ER a N tion direction is given in figure 33. It 
ra V \ is seen that above and below the plate 


/ ! \ the acceleration is downward, in front 
f LUMI \ of and behind the plate it is upward 
u | 
/ 
A 


opposite to the acceleration of the 


\ A / plate, since the air is escaping from 
ER EN, Y the plate. Lanchester asks now about 
a a ne tek NEN 7 the velocities which arise from the 

Fic. 32.— Production of circulation around a wing due to vortices leaving trailing original uniform velocity relative to 
ee the plate owing to the fact that the 


plate, while it gives rise to the accelerations as shown in figure 33, gradually comes nearer the 


air particle considered, passes by it, and finally again moves forward away from it. The pic- . 


ture of the velocities and streamlines which Lanchester obtained in this way and reproduced 
in his book was, independently of him, calculated exactly by Kutta. It is reproduced in figure 
34. It is seen that as the result of the upward accelerations of the flow away from the wing 


7 Aerodynamics I, § 110-118. 
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there is an upward velocity in front of the plate, a uniform downward acceleration at the plate 
itself due to which the upward velocity is changed into a downward one; and finally behind the 
plate a gradual decrease of the downward velocity on account of the acceleration upward. 


C. THE FINITE WING. 


16. It has been known for a long time that the aspect ratio of an aerofoil had a great effect 
on its properties. One could therefore have expected that, on account of the vanishing of 
pressure at the side edges, the intensity of the lift must decrease toward the edge, so that its 

average value for the same angle of attack must be smaller for small values 

NMW of the aspect ratio than for large ones. But the observed influence of aspect 

{ CHO + ratios is sensibly greater than could be explained in this way. We must 
ZN therefore investigate whether an explanation of this phenomenon can be found, 
Fig. 33. —Acceleration if we apply to the finite aerofoil in some proper manner the results which are 
Detay ‘lone. nat known to hold for uniplanar flow. 
plate accelerated at It is easily seen that vortices in the free fluid must here be taken into 
night angle to itssur account. For it is certain that circulation is present around the middle of 
the wing, because no lift is possible without circulation. If a closed line 
drawn around the middle of the wing, around which, therefore, there is circulation, is displaced 
sideways over the end of a wing, it will certainly no longer show circulation here when it is 
beyond the wing. From the theorem that the circulation along a closed line only changes 
if it cuts vortex filaments, and that the amount of the change of the circulation equals the 
sum of the strengths of the vortex filaments cut (see sec. 8), we must conclude that from each ` 
half of a wing vortex filaments whose strengths add up to T must proceed, which are concen- 
trated mainly near the ends of the wing. According to the Helm- 
holtz theorem we know further that every vortex produced inthe "= =—————______ 
fluid continues to move with the same fluid particles. We may — Z 
look upon the velocities produced by the wing as small compared ER ER tana 
with the flight velocity V, so that as an approximation we may — 7 LT 
assume that the vortices move away from the wing backwards 
with the rectilinear velocity V. (If it is wished, we can also im- 
prove the considerations based upon such an assumption if the 
motion of the vortices of themselves relative to the air is taken into 
account. This will, however, be seen to be unnecessary for practical applications of the theory.) 

In order now to obtain the simplest possible scheme, we shall assume that the lift is uni- 
formly distributed over the wing; then the total circulation will arise only at the ends, .and 
continue rearwards as free vortices. The velocity field of an infinitely long wing, as we saw, 

was the same at great distances as that of a rectilinear vortex 


I filament instead of the wing. We shall assume that the corre- 
| sponding statement holds for the finite wing. We thus obtain, 
<p) 


Fic. 84,—Streamlines around an infinitely 
long curved plate. 


for the velocity field around a finite wing, a picture which is 
somewhat crude, it is true, if we take for it the velocity distri- 

bution due to a vortex filament of corresponding shape. 
It may be mentioned here that, on account of there being 
A— the same laws for the velocity field of a vortex filament and the 
en magnetic field of an electric current (see sec. 8), the velocity 
ET orties tepining the ning Near a finite wing can also be investigated numerically by cal- 
culating the direction and intensity of the magnetic field pro- 


` duced near an electrical conductor shaped as shown in figure 35 due to an electrical current 


flowing in it. . 
The principles for the calculation of this velocity field have been stated in section 8; the 
total velocity is made up out of three partial velocities which are caused by the three rectilinear 


vortex portions. As is seen without difficulty, for the region between the vortices the flow is 
downward, outside it is upward. 


APPLICATIONS OF MODERN HYDRODYNAMICS TO AERONAUTICS, 


17. This approximation theorem is specially convenient if the conditions at great distances 
from the wing are treated. With its help we can explain how the weight of an airplane is 
transferred to the ground. In order to make the flow satisfy the condition that at the ground 
components of velocity normal to it are impossible, we apply a concept taken from other 
branches of physics and superimpose the condition of an image of the airplane, taking the earth 
as the mirror. On account of symmetry, then, all velocity components normal to the earth’s 
surface will vanish. If we use as our system of coordinates one attached to the airplane, we 
have then the case of stationary motion. If we take the X axis in the direction of the span 
of the wing, the Y axis horizontal in the direction of flight and the Z axis vertically down, 
and if u, v, w are the components of the additional velocity due to the vortices, then calling 
Po the undisturbed pressure and p’ the pres- 
sure difference from po, and neglecting the 
weight of the air, Bernouilli’s equation gives us 


pot p+ 5 [w+ w- V)?+-w'] = pot-5 V? 


If this is expanded and if u?, v’, and w? are 
neglected as being small of a higher order, there 
remains the simple equation 


p'=p Vv (27) 


For the determination of the pressure dis- 
tribution on the ground we must now calculate 
the value of v. Let us assume the vortices 
run off the wing in an exactly horizontal direc- 
tion (actually, their path inclines downward 
slightly), in which case they do not contribute 
tov. There remains then only the “ transverse 
vortex” of the length J (effective span) and 
the circulation T. We will assume that the 
span of the wing is small in comparison with 
the distance h of the airplane from the ground. 
In that case we can treat the transverse vortex 
as if it were a single vortex element. We ob- 
tain, then—see figure 36—at a point A, with 
the coordinates x and y, a velocity perpendicular to the plane ABF, of the amount 


Fig. 36.—Application of method of images to airplane flying near the 
ground, 


The image of the airplarfe furnishes an equal amount perpendicular to the plane ABF’ 
If Bis the angle between the plane ABF and the XY plane, then the actual velocity at the 


ground, as far as it is due to the transverse vortex, will be the resultant of v, and v,. It is there- 
1 


fore v=2ù, sin £, or, if we write sin a=, sin 8 -4 (see fig. 36) 
ja Lih (28) 
If we take into account the fact that, according to the Kutta-Joukowski formula (20), 
pT Vi=A, equations (27) and (28) lead to the relation 
-If this is integrated over the whole infinite ground surface, it is seen that the resultant 
force due to the pressures on the ground has exactly the amount A. It is thus proved that the 


23 


24 


REPORT NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS. 


pressure distribution due to the circulation motion transfers to the ground exactly the weigh 


of the airplane. The distribution of the pressure, which according to formula (29) is axially 


symmetrical with reference to the foot of the vertical line drawn from the airplane, is shown in 


figure 37. The pressure maximum is p,= — Its amount, even for low heights of flight, is 


very small, since the surface over which the pressure is distributed is very large. 

18. Applications of an entirely different kind may be made of the velocity field which belongs 
to the vortex of figure 35. For instance, an estimate may be made as to the magnitude of the 
downward velocity component at any point of the tail surfaces, and in, this manner the influence 
of the wings upon the tail surfaces may be calculated. If in accordance with the Kutta- 
Joukowski formula the lift is written A=pT VI, in which, taking account of the fact that a 


— aig Ei, 


Mee a 


FIG. 37.—Distribution of pressure on ground caused by airplane flying near it. 


portion of the vortices flow off within the ends of the wing, J, can be taken somewhat less than 
the actual span b, then at a distance d, behind the wing, the velocity component downward is 


rer i 2a- (143) (80) 


ee 

in which a= J (5) +d. 
If the flight velocity is V, this gives for the inclination of the downward sloping air-current 
tan =y We proved this relation in the year 1911 and found an approximate agreement 


with observation. 

The principle made use of above has been applied with profit to the calculation of the 
influence of one wing of a biplane upon the other wing and has given a method for the calcula- 
tion of the properties of a biplane from the properties of a single wing as found by experiments. 
The fundamental idea, which is always applied in such calculations, is that, owing to the vortex 
system of one wing, the velocity field near the wing is disturbed, and it is assumed that a wing 
experiences the same lift as in an undisturbed air stream if it cuts the streamlines of the flow 
disturbed by the other wing in the same manner as a monoplane wing cuts the straight stream- 
lines of the undisturbed flow. As is easily seen, the wing profile must in general be slightly 
turned and its. curvature slightly altered, as is shown in figures 38 and 39. By the rotation of 
the wing the direction of the resultant air force acting on it is turned through an equal angle. 
If the magnitude of the velocity as well as its direction is also changed, this must be expressed’ 
by a corresponding change in the resultant air force. 
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As an illustration we will treat briefly the case of a biplane without stagger. The most 
important component of the disturbance velocity w is again the vertical one; in the plane of 
the mean lift lines of the biplane it is affected only by the pair of vortices running off the wings, 
since the transverse vortex of one wing causes only an increase (or decrease) of. the velocity of 
flow at the other wing. . We are concerned here only with the calculation of that downward 
disturbance velocity due to the vortices from the wing not under investigation, since the other 
vortex system is present with the monoplane and its influence has already been taken into 
account in the experiments on a monoplane. 

The total velocity due to a portion of a vortex proceeding to infinity in one direction, in 
the plane perpendicular to the vortex at its end, is, as may be deduced easily from the formula 
in section 8, exactly half of the corresponding velocity in the neighborhood of a rectilinear 
vortex filament extending to infinity in both directions. This can also be easily seen from the 
fact that two vortex filaments, each extend- l 
ing to infinity in only one direction—but | 
oppositely. in.the two cases—form, if com- | 
bined, a singlefilament extending to infinity | 
in both directions. The total velocity 

caused at the point P by the vortex A, | 


LAE A 


T; 
ua = 2 2. 
see figure 40, is -> where r= YaR+ hr; its LLALL 


vertical component is 


The vertical component due to the 
vortex B is 


Thevefone the vertical component due 
to both vortices is 


-E (E+) (31) 


If we assume that the lift is uniformly 
distributed over the effective span l, which 
again we shall take as somewhat less than 
the actual span, then, since every element of 
the wing must be turned through the angle 


p according to the formula tan =p the 


direction of the air force must be turned 
also, which means a negligible change in ma , abi he otner: station o a 
the lift, but an increase in the drag of this nfluence of one wing of a bip ane upon the other; rota ion of wing profile, 
wing which must be taken into account. piteration of its curvature. 

It is essential then in this calculation that we pass from a condition for a monoplane to 
one in which the wing when part of a biplane has the same lift as when considered as a mono- 
plane. The angle of attack for which this condition will arise can be estimated afterwards 
from the average of the angles ¢. 


Tie, 39. 


19. The contribution of vortex A to the increase of the drag of the upper wing in figure. 40 
is evidently 


A 4, T nde _A, T r. 
VLO T mV e T, mVn 


W’= 
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The contribution of vortex B is, by symmetry, the same. In accordance with equation (20), 


we can put r= and thus obtain for the increase of the drag of the upper wing 
1 

_ AA, log r,/r, 

2 LL 2rp Ve 
By the symbol W,, is meant that it is the drag produced by wing 1 upon wing 2. One can 
convince himself easily that the drag W,, which wing 2 produces upon wing 1 has the same 
magnitude. Therefore the total increase of drag due to the fact that two monoplanes which 
produce the lifts A, and A, are combined to form a biplane, the two lifts remaining unchanged 

(the angles of incidence of course being changed), is 

A,A, 2 
2nl,- bg log ri 


(32) 


Wit Wy, =2 Wy= 


(33) 


in which, as always 1=3 pV?.)° 


Upon the change in the magnitude of the velocity, which in accordance with the approxi- 
mation used depends only upon the disturbance velocity v in the direction of flight, only the 
transverse vortex of the other wing has. an influence. For 

„ any point this influence, according to our formula, is given by 


v= G+") (34) 
in which r and r! have the same meaning as before. The 
> upper wing experiences due to the lower an increase in veloc- 
Z j y ity, the lower one experiences due to the upper a decrease in 
velocity, to which correspond, respectively, an increase or a 

Fic. 40.—Velocity at a point P due to the tip decrease in lift as shown by the usual formule. If we wish to 
vortices, . . . : 
keep the lifts unchanged, as required in the treatment given 
above, it is necessary to çhange. the angles of attack correspondingly. 
The effective change in the curvature? of the wing profile will, fer simplicity’s sake, be 


discussed here only for the medial plane of the biplane, i. e., for z=- It is obtained in the 
simplest manner by differentiating the angle of inclination of the air current disturbed by the 


other wing, which is, remembering that tan o=> 


dd (tan 6)= p92 (35) 


Outside of the vertical plane, owing to the disturbing wing, three vortices contribute to 
the magnitude of w. A side vortex contributes, at a point at the height A and the nuance y 


in front of the transverse plane, a velocity v” perpendicular to r’’ of the amount y5 (1-4 5, 


T 4/ 
w= -0 g 3) = =Q u a 


The transverse vortex contributes 


and therefore its share of w is 


in which the meaning of r’’ and r’’’ may be seen from figure 41. The total w is, accordingly, 


Th, /1 
w=2w, +w = IR” am = gti) 


8 The mutual action of two wings placed side by side can also be calculated from the considerations stated above, and results in a decrease of 
thedrag. This decrease is of a similar kind to that which arises in the theory of a monoplane by an increase in the aspect ratio. 

9 By change in curvature of the wing is meant that if the flow were to be kept straight and the curvature changed, the forces on the wing would 
be changed exactly as they are on the actual wing owing to the change in the fow.—Tr. 
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The differential of this with reference to y, for the value of y=o, is, since then r''=r and 


r” =h 
(20) = P(244), 
dy) — trr? ' k 


the curvature sought is, then, according to equation (35), 
1 Tr 4/1 1 
R nY P R (36) 
Calculations of the preceding nature were made in 1912 by my assistant, A. Betz, so as 
to compare experiments with monoplanes and biplanes and to study the influence of different 
angles of attack and different degrees of stagger of the two wings of a biplane upon each other. 
The influence upon the drag was not known to us at that time, and the calculation was carried 
out so as to obtain the changes in the lift due to w, to v and to the curvature of the streamlines. 
In this connection the change of the lift of a monoplane when flying near the earth’s surface 
was also deduced, by calculating the influence of the “mirrored wing” exactly as was that of 
the other wing of a biplane. All that was necessary was to change some algebraic signs, 
because the mirrored wing had negative lift. The theory. of these 
calculations was given by Betz in the Z. F. M., 1914, page 253. 
The results of the theory of Betz, from a more modern 
- standpoint, such as adopted here, were given in the Tech- 
nische Berichte, volume 1, page 103 et seq. There one can 
find the discussion requisite for the treatment of the most 
general case of a biplane having different spans of the two 
wings and with any stagger. In the case of great stagger 


. . . $ Fig. 41.—Curvature of wing-profile at 
it appears, for example, that the forward wing is in an ascend- its middle point due to velocities 


ing air current caused by the rear wing; the latter is in an used by transverse and tip vortices. 
intensified descending current due to the forward wing and the vortices flowing off from it. 
Corresponding to this, if the angle of attack is unchanged, the lift of the forward wing is 


increased, and that of the rear one weakened; at the same time the ratio Drag experiences & 


decrease for the forward wing and a marked increase for the rear one. 

For a wing in the neighborhood of the ground, owing to the influence of v there is a decrease 
of lift, and conversely there is an increase of lift due to the influence of w, provided the angle 
of attack is kept constant, but as the result an evident decrease in the ratio Drog, Owing to 
this last it is seen why in the early days of aeronautics many machines could fly only near the 
ground and could not rise far from it. Their low-powered engines were strong enough to over- 
come the diminished drag near the ground but not that in free air. 


D. TEEORY OF THE MONOPLANE, 


20. If we extend the principles, which up to this point have been applied to the influence 
of one wing upon another, to the effect upon a single wing of its own vortices, it can be said in 
advance that one would expect to find in that case effects similar to those shown in the influence 
of one wing of a biplane upon the other, i. e., the existence of lift presupposes a. descending 
flow in the neighborhood of the wing, owing to which the angle of attack is made greater and 
the drag is increased, both the more so the closer to the middle the vortices flowing off at 
the ends are, i. e., the smaller the aspect ratio is. One might propose to apply the theory pre- 
viously given for biplanes by making in the formulas of this theory the gap equal to zero. Apart 
from the fact that the formulas developed do not hold for the immediate neighborhood of the 
vortex-producing wing, but must be replaced by more accurate ones, this certainly is not the 
proper path to follow, for, in the earlier treatment, we have taken, the undisturbed monoplane 
as the object with which other cases are to be compared and have asked what drag, what change 
in angle of attack, etc., are caused by adding a second wing to this monoplane. .To proceed 
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according to the same method, we must seek for the theory of monoplanes another suitable 
object of comparison. As such, the infinitely long wing will serve. Where the discussion 
previously was about change of angle of attack, increase of drag, etc., we intend now to refer 
these to the infinitely long wing as a starting point. Since in the theoretical nonviscous flow 
the infinitely long wing experiences no drag, the total drag of such a wing in such a fluid must 
be due to vortices amenable to our calculations, as the following treatment will show. In a 
viscous fluid drag will arise for both wings; infinitely long or not, which for those angles of 
attack for which the profile is said to be ‘‘good” is, according to the results of experiment, of the 
ordgr of magnitude of the frictional resistance of.a plane surface. 

The carrying out of this problem is accompanied with greater difficulties than the calcula- 
tion for a biplane as given. In order to obtain the necessary assistance for the solution of the 
problem, we shall first be obliged to improve the accuracy of our picture of the vortex system. 

The density of the lift (lift per unit length) is not constant over the whole span, but in 
general falls off gradually from a maximum at the middle nearly to zero at the ends. In ac- 
cordance with what has been proved, there. corresponds to this a circulation decreasing from 
within outward. Therefore, according to the theorem that by the displacement of the closed 
curve the circulation T can change only if a corresponding quantity of vortex filaments are cut, 
we must assume that vortex filaments proceed off from the trailing edge wherever T changes. 


For a portion of this edge of length dx the vortex strength is therefore to be written ar dz, and 
hence per unit length of the 
edge is a These vortex fila- 


. ments flowing off, closely side 
by side, form, taken as a whole, 
a surface-like figure, which we 
shall call a ‘‘ vortex ribbon.” 

For an understanding of 
this vortex ribbon we can also 
‘approach the subject from an 

entirely different side. Let us consider the flow in the immediate neighborhood of the surface of 
the wing. Since the excess in pressure below the wing and the depression above it must vanish 
as one goes beyond the side edges of the wing in any manner, there must be a fall in pressure near 
these edges, which is directed outward on the lower side of the wing and inward on the upper. 
The oncoming flow, under the action of this pressure drop, while it passes along the wing, will 
receive on the lower side an additional component outward, on the upper side, one inward, 
which does not vanish later. If we assume that at the trailing edge the flow is completely 
closed again, as is the case in nonviscous flow, we will therefore have a difference in direction 
between the upper and lower flow; the upper one has a relative velocity inward with reference 
to the lower one, and this is perpendicular to the mean velocity, since on account of the Ber- 
nouilli equation in the absence of a pressure difference between the two layers the numerical 
values of their velocities must be the same. This relative velocity of the two flows is exactly 
the result of the surface distribution of vortices mentioned above (as the vortex theory proves, 

a surface distribution of vortices always means a discontinuity of velocity between the regions 

lying on the two sides of the surface). The relative velocity is the greater, the greater the side- 

wise pressure drop, i. e., the greater the sidewise change in lift. The picture thus obtained 
agrees in all respects with the former one. 

21. The strengths of our vortex ribbon remain unchanged during the whole flight, yet. 
the separate parts of the ribbon influence each other; and there takes place, somewhat as is 
shown in figure 42, a gradual rolling up of the ribbon, as a closer examination proves. An 
exact theoretical investigation of this phenomenon is not possible at this time; it can only be 
said that the two halves of the vortex ribbon become concentrated more and more, and that 
finally at great distances from the wing there remain a pair of vortices with rather weak cores. 


Fıe. 42.—-Change in shape of vortex ribbons at great distances behind the wing. 
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For the practical problem, which chiefly concerns us, namely, to study the reaction of the 
vortices upon the wing, it is not necessary to know these changes going on at a great distance, 
for the parts of the vortex system nearest the wing will exercise the greatest influence. We 
shall therefore not consider the gradual transformation of the vortex ribbon, and, in order to 
make the matter quite simple, we shall make the calculation as if all the vortex filaments were 
running off behind in straight lines opposite to the direction of flight. It will be seen that, 
with this assumption, the calculations may be carried out and that they furnish a theory of 
the monoplane which is very useful and capable of giving assistance in various ways. 

If we wish to establish the method referred to with greater mathematical rigor, we can 
proceed as follows: Since the complete problem is to be developed taking into account all cir- 
cumstances, we shall limit ourselves to the case of a very small lift and shall systematically 
carry through all calculations in such a manner that only the lowest power of the lift is retained, 
all higher powers being neglected. The motion of the vortex ribbon itself is proportional to 
the total circulation, therefore also proportional to the lift; it is therefore small if the lift is 
small. If the velocities caused by the vortex ribbon are calculated, first for the ribbon in its 
actual form, then for the ribbon simplified in the manner mentioned, the difference for the two 
distributions will be small compared with the values of the velocity, therefore small of the 
second order, i. e., small as the square of the circulation. We shall therefore neglect the differ- 
ence. Considerations of this kind are capable of deciding in every case what actions should be 


taken into account and what ones may be neglected. © By our simplifications we have therefore. 


made the problem linear, as a mathematician says, and by this fact we have made its solution 
possible. It must be considered a specially fortunate circumstance that, even with the greatest 
values of the lift that actually occur with the usual aspect ratios, the independent motion of the 
vortex ribbon is still fairly small, so that, in the sense of this theory, all lifts which are met in 
practice may still be regarded as small. For surfaces having large chords, as, for instance, a 
square, this no longer holds. In this case there are, in addition, other reasons which prove 
that our theory is no longer sufficiently accurate. This will be shown in the next paragraph. 

It has already been mentioned that the infinitely long wing will serve as an object of com- 
parison for the theory of the monoplane. We shall formulate this now more exactly by saying: 
Every separate section of the wing of length dz shall bear the same relation to the modified flow 
due to the vortex system as does a corresponding element of an infinitely long wing to the recti- 
linear flow. The additional velocities caused by the vortex system vary from place to place 
and also vary in the direction of the chord of the wing, so that again we have to do with an 
influence of curvature. This influence is in practice not very great and will for the sake of 


simplicity be neglected. This is specially allowable with wings whose chords are small in. 


comparison with their spans, i. e., with those of large aspect ratio. If one wishes to express 
with mathematical exactness this simplifying assumption, it can be said that the theory of an 
actual wing of finite chord is not developed, but rather that of a “lifting line.” It is clear 
that a wing of aspect ratio 1:6'may be approximated by a lifting line, specially if one considers 
that actually the lift is concentrated for the greatest part in a region nearer the leading edge. 
It is easily seen, however, that a surface in the form of a square can be approximated only 
poorly by a lifting line. 

If we assume a straight lifting line, which lies in a plane perpendicular to the direction 
of flight, the flow due to the vortices, which according to the Biot-Sarvart law, is caused by its 
own elements, will not produce any velocities at the lifting line itself except the circulation 
flow around it, which would also be present for an infinitely long lifting line having the same 
circulation as at the point observed. All disturbance velocities at a point of the lifting line, 
which are to be looked upon as deviations from the infinitely long lifting line, are due therefore 
to the vortices which run off and hence can be calculated easily by an integration. 

A qualitative consideration of the distribution to be expected for the disturbance velocities 
along our lifting line shows at once that—just as was the case for a biplane—the chief thing is the 
production of a descending current of air by the vortices. If we wish to retain the lift of the 
same intensity as with the infinite wing, the angle of attack must be increased, since the descend- 
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ing air stream added to the wind due to flight causes 


a velocity obliquely downward. In addi- 


tion, the air force, as before, must be turned through the same angle, so that a drag results. 


The rotation will be the greater, the greater the lift 


w 


/ 


Fig, 43.— Wing having finite, but small, chord. Distribution of vertical velocity 
component along a line parallel to direction of flight. 
I. Infinitely long wing. 
IT. The downward velocity produced by vortices flowing off. 
III. Finite wing, sum of I and IT. 


and the closer to the middle of the wing 
the main production of vortices is. The 
drag must therefore increase both with 
increasing lift and with decreasing span. 
A picture of what occurs with a 
wing of finite but small chord is given in 
figure 43. There the change is shown 
of the vertical velocity component along 
a straight line parallel to the direction 
of flight through the middle of the 
wing; in the upper part of the diagram, 
for the infinitely long wing, in the lower 
part, for the finite wing. We see from 
Curve I the rising flow in front of the 
wing, its transformation into a descend- 
ing one at the wing itself and the gradual 
damping of the descending component 
due to the upward pressure drop behind 
the wing. (See sec. 15.) The corre- 
sponding curve for the finite wing is 
Curve III. It is derived from I by 
adding to the latter the descending ve- 
locity II. We recognize the rotation of 
the profile as well as that of the lifting 
force, which was originally perpendicular, 


through the angle ¢ where tan ¢ -7 and 


w is the velocity downward at the location of the center of pressure (i. e., at the lifting line). 
If we follow the method of Lanchester, as described in section 15, the downward velocity w can 
also be looked upon as a diminution of the ascending flow at the leading edge of the wind due 
to the absence of the sidewise prolongation of the wing, i. e., to the deviation from an infinitely 
long wing which was the basis of the treatment in section 15. Dis- 


cussions very similar to this are given in Lanchester, Volume I, 


section 117. 


It may be seen from the figure that at great distances behind 
the wing the descending velocity is 2w, which agrees with the relation 
already mentioned that the velocities due to a straight vortex fila- 
ment extending to infinity in both directions are twice those due to 
a filament extending to infinity in one direction only, for points in the 
plane perpendicular to this vortex passing through its end point. 

22. The mathematical processes involved in carrying out the 
theory outlined above become the most simplified if one considers as 


known the law, according to which the lift is distributed over the wing. Fie- 44-Velocity at a point ~ due 


We shall call this the “first problem.” The calculation is made as 


to vortex leaving wing at point x. 


follows: The distribution of lift is the circulation expressed as a function of the abscissa x. The 


strength of the vortex filament leaving an infinitely small section dz is then a - de. This 


roduces at a point x’, according to what has been already explained, a vertical velocit 
p P g y exp ’ y 


downward or upward of the amount J 
1 gr 


dx 


= a’ wa 
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In this z’—2 takes the place of r in section 8. If the circulation falls to zero at the ends of the 
wing, as is actually the case, then all the vortices leaving the wing are of this kind. The whole 
added velocity at the position 2’, rn that the function T(x) is everywhere continuous, is 


ddr ie 


o de rg (87) 


ay 


We must take the so-called “chief eine of the integral, which is indeterminate at the 


point =g’, i. e., the limiting value a ‘ 
lim. ( 1 2 | ) 
e—0 o THe 


must be formed, as a closer examination shows. We can do this by calculating, instead of the 
value of the velocity at the lifting line, which is determined by the preponderating influence 
of the nearest elements, the value of w for a point a little above or below the lifting line. It 
is seen that this last is not indeterminate and that by passing to a zero distance from the lifting 
line it reaches the above limit. Concerning this excursus, important in itself, the preceding 
brief remarks may be sufficient. 

After the calculation of the integral of (37), the downward velocity is known as a function 
of the abscissa x’ (which > later shall call x). We then also know the inelination of the 


resultant air flow, tan ¢= y the lift dA = T dx’, acting on the section dg’, therefore con- 
tributes to the value of the drag 

dW=tan p - dA =p (2) - wir’) de’ 
since it is inclined backward by the small angle y. The total es is therefore 


Ten? a 
Wap | r(')-w(e") de! = |) zen a 


For a long time it was difficult to find suitable functions to express the distribution of lift, 
from which a plausible distribution of w would be obtained by equation (37). After various 
attempts it was found that a distribution of lift over the span according to a half ellipse gave 
the desired solution. According to this, if the nn of coordinates is taken at the center of the 


wing, Tr 
a in OR — 
r=. EEA (v5 D)o hence 4 = Se 


ay(a)-# 
The “chief value” of the integral 


tdt 
CE is equal to + 
and therefore the integral of equation (37) is equal to py and thus is independent of 2’ 
and constant over the whole span. Hence 


WOR 
The value of T, is obtained from 


a +b/2 a +b/2 ae i , E 7 
vf Pde oyn f V- de =pV V4 b, 


giving 
pa tA 
° aepVb 
Hence 
2A 
= ToVE (39) 
Since w is constant there is no need of calculating the drag by an integral, for it is simply 
w 2A? A? 
Weed (40) 
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The calculation can also be performed for distributions of circulation given by the following 
general formula: 


T=Vi-F (Mot TPT, B+ :.) (41) 
in which £= TA l 
According to the calculations of A. Betz 
m=n 
ES 2 [(2n + 1) Pa-m— 2NPn-m 4 (42) 
z m=o 
and Me 
W= "2, Zy | Tzi Dox 2, Gism [(2k +1) Pr_m— Bin (43) 
m=o 
in which the numbers p and q have the meaning 
1.3. en = 
Pa= TF. > a= or Po=l, p50 


The elliptical distribution of litt, a from its simplicity, has obtained a special meaning 
from the fact that the drag as calculated from equation (40) proved to be the smallest drag that 
_ is imaginable for a monoplane having given values of 
the total lift, the span and the velocity. The proof of 
this will be given later. 

It was desirable to compare this theoretical mini- 
mum drag with the drags actually obtained. As far back 
as 1913 this was done, but, on account of the poor quality 
of the profiles then investigated, all that was done was 
to establish that the actual drag was greater than the 
theoretical. Later (1915) it was shown, upon the in- 
vestigation of good profiles, that the theoretical drag 
corresponds very closely to the relation giving the change | 
of the observed drag as a function of the lift. If we 
plot in the usual manner the theoretical drag, as given 
in formula (40) as a function of the corresponding lift, 
we obtain a parabola, which runs parallel with the 
measured “polar curve” through the entire region for 
which the profile is good. (See fig. 45.) 

This process was repeated for wings of different as- 
pect ratios, and it was proved that for one and the same 
profile the difference between the measured and the theo- 

Fig. 45.—Polar diagram showing theoretical drag and retical drags for one and the same value of the lift coeffi- 
ER cient had almost identically the same value in all cases. 

This part of the drag depends, however, upon the shape of the profile, and we have therefore 
called it “profile drag.” The part of the drag obtained from theory is called “edge drag,” 
since it depends upon the phenomena at the edges of the wings. More justifiably the expres- 
sion “induced drag” is used, since in fact the phenomena with the wings are to a high degree 
analogous to the induction phenomena observed with electric conductors. 

Owing to this fact that the profile drag is independent of the aspect ratio, it became possible. 
from a knowledge of the actual drag for one aspect ratio to calculate it for another. To do 
this, we pass from the formula (40) for the drag to the dimensionless lift and drag coefficients, 


~ 
by letting i =¢, and mn we obtain then for the coefficient of the induced drag the 


relation 
2F 
[tz (44) 
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Fig, 46,— Polar diagrams for seven wings, aspect ratios 1:7, Tig. 47.--Lift coefficients plotted as function of angle of attack for aspect 
1:8, ete., 1:1. ratios 1:7 to 1:1. 
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Fic. 48.—Polar diagrams reduced from observations on Fig. 49.—Lift coefficients as function of angle of attack, reduced for 
aspect ratio 1:5. aspect ratio 1:5. 
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The profile drag may then be written ¢yo=¢y—¢wi- If this drag coefficient depends only 
upon the lift coefficient, then it would be evident, since it would be the difference between the 
measured and the theoretical drag coefficients, that, for the polar curves of two different wings 
having Ca = Cay = Ca; 

e EE e Cale 
wi rbe OW? ab” 
and therefore 


c@/F, F, 
Cw = Cw yt T (53 p2 (45) 


In a similar manner a calculation for the angle of attack may be made if we presuppose an 
elliptical distribution of lift. According to our assumption there is a close connection between 
the lift of the separate elements of the wing and the “effective” angle of attack, which is the 
same as the angle of attack of an infinitely long wing. This effective angle of attack, according 
to our earlier considerations, is the angle of attack of the chord with reference to the resultant 
air current. Itis therefore a’=a—¢. If we substitute tan & =y for ¢, and introduce in equa- 
tion (39) again the lift coefficient, instead of using the lift, we obtain for the comparison of two 

wings, expressing the fact that the effective angle of attack 
5 to be the same for two equal lift coefficients, the relation 


E A à F 


a =a. 2 
1 6? 2 Tr b? 


which leads to the transformation formula 


l . Al: + (4: (46) 
Fia. 49ya.—Uniplanar flow in case of elliptic dis- 2 1 

tribution of lift on a lifting line. These formulas have been found to hold for distribu- 
tions of lift which do not deviate too much from elliptical ones, although strictly speaking they 
apply only to the latter. The fact that the type of distribution does not have a marked effect 
is based upon the consideration that both in the calculation of drag and in that of the mean 
effective angles of attack we are concerned with average results. For the calculation of the drag 
one can also introduce the thought that no quantity varies much in the neighborhood of its 
minimum. Closer investigation of the square cornered wing has shown that, if the aspect ratio 
is not too small, the lift distribution does not deviate greatly from the elliptic type, and that the 
theoretical drag for usual aspect ratios at the most is 5 per cent greater than for the elliptic 
distribution. As an example of these formulas we shall take ‘our figures from the book published 
by the Géttingen Institute (Ergebnisse der Aerodynamischen Versuchsanstalt, 1, Lieferung, 1921). 


The first and second figures show the polar curves, and the connection between lift coefficient 


and angle of incidence for seven wings of aspect ratio’ 1:7 to1:1. The last two figures give the 
results of calculating these experimental quantities from the results for the wing having the aspect 
ratio 1:5. It is seen that, apart from the aspect ratios 1:1 and 1:2 practically no deviations 
are present. The fact that the square can not be correctly deduced from the aspect ratio 1:5 
need not excite surprise, since the theory was developed from the concept of the lifting line, and 
a square or a wing of aspect ratio 1:2 can scarcely be properly approximated by a lifting line. 
On the other hand it is a matter of surprise that an aspect ratio of 1:3 can be sufficiently 
approximated by the imaginary construction of a lifting line. The deviations in the case of the 
square are moreover in the direction one would expect from a lift distribution expanded over 
the chord. A quantitative theory is not available in any case at the present time. 

28. If the lift distribution is not given, but, for.example, the downward velocity, then the 
method of treatment followed hitherto may be used, by developing the downward velocity in a 
power series and determining the constants of the series given above for the lift from the constants 


10 The American practice is to define aspect ratio as the ratio of span to chord, which would involve taking the reciprocals of the ratios given 
inthe text. Tr. 
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of this power series, by the solution of linear equations. By this the lift distribution and every- 
thing else are known. 

“Another method for the solution of this “second problem” will be obtained by the fol- 
lowing consideration: The velocities at a distance behind the wing, on account of the connection 
mentioned so often between a vortex filament extending to infinity in one direction only and 


one extending to infinity in both directions, are twice as great as those in the cross section of. 


the lifting line, if we do not take into account the change in shape of the vortex ribbon. We 
therefore have here, neglecting this change in shape, an illustration of a two-dimensional fluid 
flow (uniplanar flow), for which the vertical velocity components at the point where the wing 
is reached are specified. For the simple case that the vertical velocity w is constant, as was 
found to be true for the elliptical lift distribution, the shape of the flow that arises has been known 
for along time. Itis given in figure 49a. It is the same as that already considered, in another 
connection, in section 15. The picture of the streamlines show clearly the velocity disconti- 
nuity between the upper and lower sides of the vortex ribbon, indicated by the nick in the stream- 
lines, and also the vortical motion around the two extreme points of the vortex ribbon, corre- 
sponding to the ends of the wing. 

Any problems of this kind can therefore be solved by means of the methods provided by 
the potential theory for the corresponding problem of two-dimensional fluid flow. We can 
not go into these matters more closely at this time; by a later opportunity some special rela- 
tions will be discussed, however. 

A “third problem” consists in determining the lift distribution for a-definite wing having a 
given shape and given angle of attack. This problem, as may be imagined, was the first we 
proposed; its solution has taken the longest, since it leads to an integral which is awkward to 
handle. Dr. Betz in 1919 succeeded after very great efforts in solving it for the case of a square- 
cornered wing having everywhere the same profile and the same angle of attack. The way the 
solution was obtained may be indicated briefly here. We start, as before, from the relation 


=g FREE EIN 4 2 
a= a + Oa! +47 


By equation (37) w is expressed in terms of the circulation. The effective angle of attack 
a’ can be expressed in terms of T, since, according to the assumptions made before the lift, 
distribution, which is proportional to T, depends directly upon a’. The relation between a’ 
and T can be assumed to be given sufficiently exactly for our purposes by a. linear expression 


T= Vt (ce, a’ +¢,) (47) 


in which ¢ is the length of the chord (measured in the direction of flight). By the introduction 
of the factor Vt, c, and c, are made pure numbers. The numerical value of c,, which is the more 
important, can be expressed, if ca, is the lift coefficient for the infinitely long wing at the angle 
of attack «’, by the relation 


1 Wary 
der 
In fact A yI 
Caw = ae 2 en (cya! +c) 
I it pV? 


For a flat-plate theory proves that c,=r, for curved wings it has a slightly greater value. 
If, according to what has gone before, we express a’ by T and w by A and, write 
= =f(x) and therefore T = ING da 
we obtain after a simple calculation the integral equation 


f k fader É [lee = Vt (ca +¢,) = const. (48) 
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A solution of this equation can be obtained by expanding I as in equation (41) and develop- 
ing then all the expressions in power series of A For every power of é there is then 


a linear equation between the quantities To, T,, etc. There is a system of equations, then, 
with an indefinite number of equations for an infinite number of unknowns; the solution of 
which in this form is not yet possible. The aspect ratio of the wing appears in these equations 


as a parameter, and it is clear that the solution for a small aspect ratio, i. e., u is easier than for 


a large one. Dr. Betz proved that a development in powers can be’made for the unknowns in 
terms of a parameter containing the aspect ratio. The calculations which are contained in the 
dissertation! of Dr.. Betz (1919) are 
very complicated and can not be re- 
produced here; but certain results will 
be mentioned. The Betz parameter 
L has the meaning 


In the application to surfaces which 
are investigated in wind tunnels the 


‘ de, . dla op 
value da E known, not re For 


this case theory gives a relation which 
can be expressed approximately 
b da 


L=3.85 T de, 1.3. 
, 


We can thus obtain the value of a 


an 


from the connection mentioned. 
| The distribution of lift density 
ag a 02 09 04,05 06 07 08 03 10 ever the span iselliptical for verysmall 
Fic. 30.—Change in distribution of lift, as a function of L, the parameter of Betz. aspect ratios and for greater ratios 
becomes more and more uniform; for 
very elongated wings it approaches gradually a rectangular distribution. Figure 50 shows 

this change in the distribution depending upon L. 

The drag of the wing with the rectangular distribution is greater naturally than with the 
elliptic distribution, since this gives the minimum of drag, yet the differences are not very great; 


for instance, for L=4, (;=6) it is about 5 per cent greater than that of the elliptical distri- 


bution. An approximation formula, according to the values obtained by Betz, is 
OR 
We zB (0.99 +0.015 L). 


This is applicable for values of L between 1 and 10. 

The distribution of lift, downward velocity, and drag upon a very elongated wing is shown 
qualitatively in figure 51.. It is seen that the downward velocity and the drag gradually accumu- 
late at the ends of the wing. This gives also the correct transition to an infinitely long wing, 
with which for interior positions the lift is constant, and the downward velocity and the drag 
are equal to zero, while, as we know, near the ends these last quantities always assume finite 
values. 


u Printed in extracts in Heft 2 of the “ Berichte u. Abhl. der Wiss. Ges. f. Luftf.” Munich, 1920 (R. Oldenbourg). 
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E. IMPROVED THEORY OF AIRPLANES HAVING MORE THAN ONE WING. 


24. The knowledge obtained in the theory of a monoplane can be applied also to multi- 
planes and furnishes here a series of remarkable theorems. We shall limit ourselves to the 
theory of the first order, as designated in the theory of monoplanes, therefore we shall neglect 


the influences of v. Further, we shall not take into account the effect of curvature—i. e., we - 


shall consider the separate wings replaced by “lifting lines.” For the sake of simplicity we 
shall limit ourselves to multiplanes with wings which are straight and parallel to each other. 
The generalization of the theorems for nonparallel wings, corresponding to the deduction given in 
“ Wing theory II,” will then be stated without proof. 

Let us first solve the introductory problem of 
calculating the vertical velocity w produced by a 
lifting line at a point A which lies off the lifting 
line. At the beginning let us assume that this point 
lies in the same “ transverse plane” (plane perpen- 
dicular to the direction of flight). According to 
our assumption as to the location of A, the action 
of the transverse vortex is zero: With reference to 
the longitudinal vortices it is to be remembered that 


| | 
the velocity q- > & produced by a longitudinal m (mmm 


OD fg av p Pe te F 
vortex of strength g T is perpendicular to the Fie. 51.—Distribution oflift, down-wash and drag for along wing. 


line a (see fig. 52), and therefore must be multiplied by sin 8 to obtain the vertical component. 
We arrive at the downward velocity, therefore, by integrating over the lifting line, viz: 


b 
_ 1 fdT /sin ß 
ve) = 


This relation can be brought into another form by a partial integration. Since at both wing 


ends TF =0, we have 
_1 (>, d /sin u 
vr {r dz ( a 7 


But 
d (sin B)_ @ = «—2¢? [1—2 sin? B_ cos 28 
: d\ a ) dt\@) a ~ a? a 
so that we have 
_1 f®,cos28 
Ber dz (50) 


With the aid of this relation we can write down immediately the value of the drag which arises 
owing to a second wing being under the influence of the disturbance caused by the first wing 
lying in the same transverse plane. Let us call w,, the disturbance velocity at 


A #__ a point A-on the second wing. According to the results of section 22, the 
drag then is 
þ2 
K S W= p [Twade 
© -= ` or, if the value of w, as given by equation (50) is substituted, 
h 52 i 8 he bi bs cos 26 1) 
peak at the Wand fy SE dade a 


lifting line, but in . 5 . å iya 
the transverse The double integral, as one sees, is perfectly symmetrical in the quantities asso- 


plane, due to the ciated with both wings 1 and 2. We conclude from this that the drag which 
ae wing `l experiences owing to the presence of wing 2 is of the same amount as 
the drag calculated here, that, therefore, 
Wa= Wa: 
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In the more general case of two curved lifting lines lying in a transverse plane, a formula is 
obtained which differs from equation (51) only in having cos (8, +£) in place of cos 2 8 in which 
8, and £, are the angles which the line a makes with the normals on the two lifting elements con- 
nected by a, and in having ds,ds, in place of dz,dz,. The relation W,,= W,, therefore holds in 
this case also. This mutual relation, which was discovered in a different manner by my assist- 
ant, Dr. Munk, is of importance in various applications. Since it plainly is not necessary for 


the lifting elements taken as a whole to belong to a single surface, the theorem may be stated: 


If, out of a lifting system all of whose elements lie in a transverse plane any two groups 
are selected, the portion of the drag experienced by group 1 due to the velocity field of group 2 
is exactly of the same amount as that experienced by group 2 due to the velocity field of group 1. 
We can interpret the partial integration performed above by saying that the velocity w 
appears by it as built up out of the contributions by merely infinitesimal wings having the 
length dz and the circulation T, while previously we have always built it up out of the actions 


of the separate vortices Gide, The integrand of equation (50) in fact agrees with the velocity 


which is caused by two vortex lines of equal but opposite strengths T lying at a distance dx apart. 
The double integral in equation (51) can, from this point of view, bé looked upon as the sum 
of the actions of the vortex strips of all the elements dz, on all the 
lifting elements dz,. 

The objection might be raised that equations (50) and (51) cease 
to be applicable if the value a=o appears, since this gives an expres- 
sion of the form o-—0o. They are not, therefore, suited for the 
calculation of the velocity w at the wing itself. In this case we must 
return to equation (49), and take the “chief value” of the integral; 
or, the value of w, and of W, for a lifting line that lies very closc 
must be calculated, and then we can obtain our final result by passing 

ass diate eee: OO the limit for coinciding lifting lines. As is seen from this, the 
not inthetransverseplane,due Telations W,,= W,, hold also for lifting lines ‘coinciding in space, 
to the vortex system. which, besides, may have any arbitrary lift distribution. 

The mutual drag need not, as has already been mentioned, always be positive. For in- 
stance, it is negative for two wings placed side by side, since then each wing is in an ascending 
current caused by the other, and the total drag is therefore less than the sum of the mutual 
drags which each of the wings would have at a greater distance apart. The behavior of certain 
birds which in a common flight space themselves in a regular phalanx can be explained by 
reference to this. 

25. In order to be able to treat the case of staggered wing systems, the next problem is 
to calculate the velocity field due to a lifting element of the length dx together with its pair of 
vortices at a point A which may now lie off the transverse plane, and at a distance y from it. 
(See fig. 53.) The origin of coordinates will be taken at the projection of the point A upon 
the transverse plane, and the X axis parallel to the direction of the element. ` Using the abbre- 
viations 


@=242, r= +y, 
the velocity produced at the point A by one of the two vortices, by formula (6b), is given by 


adit! A the component in the direction of the Z axis, to which we here again limit our- 


u-2,(1 +2) 


The pair of vortices produces then a velocity which may be written as the difference of the 
effects of two vortices which are close together: 


Ow, Tdz[ (a? — 22? +% 
Te ea | at (1+ a = | 


selves, is, then, putting sin B=% 
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To this must be added the contribution of the transverse vortex 


Tax. 
Er Do 
The sum of these two velocities, if the two angles defined in figure 53 are introduced, amounts to 
Tdx { cos 28 : _ Sina . cos? 
dw =| a (1+sin a) + Sn se (52) 


With the help of this formula we can now calculate at once the drag experienced by a lifting 
element situated at the point A and parallel to the former, whose length is dz, and circulation 
T,. If the first element is given the index 1, this drag is 
] 2 
p Wy eats [ {Ss 28 (+ sin a) +52 a cos £] ; (53) 


72 


As is easily seen, the drag produced on the lifting element 1 by the lifting element 2 is obtained 
if in place of a and £8 the values a +r and 8 +r are introduced. Therefore it is 
2 

g Wp = 2T dar de, es, EEE sin a gos% ] (3a) 
It is seen from this that the two a of the drag are equal only if a=o, that is if the two ele- 
ments lie in the same transverse plane. Yet in the.general case the sum d? W,,+d? W,, is inde- 
pendent of a, therefore independent of the amount of stagger. The sum of the two mutual 
drags leads thus to the same formula as that already derived for nonstaggered wings. If we 
again pass to the general case of nonparallel lifting lines, in which again ds, and ds, are to be 
written in place of dr, and dz, we obtain as may be proved by performing the calculation, the 


relation 
D ` T, ds, d N B 8) 
e $3 i Qn Sf = i g l ; x (54; 


As is evident, this sum remains unchanged if the two lifting groups are displaced in the direc- 
tion of flight. Since the total drag of a lifting system is composed of such mutual drags as cal- 
culated above and of the proper drags of the separate wings, which likewise are not changed 
by a displacement of the wing in the direction of flight, the following theorem may be stated: 

The total drag of any lifting system remains unchanged if the lifting elements are displaced 
in the direction of flight without changing their lift forces. 

This “stagger theorem” was likewise proved by Munk. For a proper understanding of 
this theorem it must be mentioned expressly that, in the displacement of the separate lifting 
elements, their angles of attack must so be changed that the effective angles of attack and there- 
fore the lifting forces remain unaltered. 

This theorem, which at first sight is surprising, may also be proved from considerations of 
energy. Let us remember that, by the overcoming of the drag, work is done, and that in a non- 
. viscous fluid, such as we everywhere assume, this work can not vanish. Its equivalent is, in 
fact, the kinetic energy that remains behind in the vortex motions in the rear of the lifting 
system.. This energy depends only upon the character of these vortices, not upon the way in 
which they are produced. If we neglect, as we have throughout, any change in shape of the 
vortex system, then, in fact, the staggering of the separate parts of the lifting system can not 
have any influence upon the total drag. 

26. For the practical calculation of the total drag of a multiplane, we have then the follow- 
ing: The total drag consists of the sum of all the separate drags and of as many mutual drags 
as there are combinations of the wings in twos. If the nature of the lift distribution over all 
the separate wings is specified, then the proper drags are proportional to the square of the sepa- 
rate lifts; the mutual drags, to the product of the lifts of the two wings in question. If the 
coefficients of this mixed quadratic expression are all known, then one can solve without diffi- 
culty the problem: For a specified total lift, to determine the distribution of lift over the sepa- 
rate wings which will make the total drag a minimum. 
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In order to know these coefficients to a certain degree I calculated them for the case of two 
straight lifting lines whose middle points lie in the same plane of symmetry, with the assump- 
tion that the lift over each separate wing is distributed according to a half ellipse. The results 
are given in my paper “The induced drag of multiplanes” in Volume III, part 7 of the Tech- 
nische Berichte. For this purpose the velocity w for the entire neighborhood of a wing in the 
transverse plane was first calculated by formula (49), and then the integrals for the mutual 
drags were obtained by planimetry. To show the analogy with equation (40) this may now 
be expressed by the formula eu 

es fag h (55) 


by means of which the numerical factor e can be expressed as a function of the two variables 


-a and = Calculation gave the following table: 
1 


>, + b,) 
TABLE I.— Values of o 
| | | 
| 2h 0.05 0.1 | 015 0.2 0.3 0.4 0.5 
bitb: | à i 

by [18 0.780 | 0.655 | 0.561 | 0.485 | 0.370 | 0.290 | 0.230 

br) lg soo | :600 | 523 | 2489 | 7355 | 282 | 225 

1 | 26 1540 | 1485 | 1437 | 1394 | 7315 | 1255 | 1210 


The curve of the function ø is given in figure 54. For the most important case, viz, for 
two wings of equal span, I have developed an approximation formula which is 


1— 0.667 
a 5 (56) 
1.055 +3.77 


It may be used from R 0.05 to 0.5. 


The total induced drag of a biplane is then, if b, is the greater span and if the ratio = 
is designated by u : 


W= Wi +2 Wiz + Wy en = 20uA,A, F pA) (57) 


Simple calculation shows that for a given A; +4, this drag is a minimum for 


1 

A, : Ay= (4-0) :(2-0) (58) 

The value of the minimum is fouud to be 

_(A,+4,)? 1-0? 

Wanin = ngb 1-Qepte (59) 
The first factor of this formula is the drag of a monoplane having the span b, and the lift A, +4,. 
Since e< p, the second factor of the formula is always less than 1, i.. e., the induced drag of a 
biplane is less than that of a monoplane which in the same span carries the same load. For 


‘ a “tandem,” i. e., an arrangement of two wings one behind the other, the stagger theorem shows 


an equivalence with two coinciding wings, i. e., a monoplane. Among the different biplanes 
having prescribed span b, and prescribed gap h, that one is the most favorable in which the 
second wing also has the span b,. The most favorable ratio of the two lifts is then 1:1 and 


the second factor of equation (59) becomes equal to 5 (1+0). 


These statements must not, however, be misunderstood; they refer only to the comparison 
of such wing systems as have the same value for the greatest span. Naturally, for every biplane 
a monoplane may be found with somewhat greater span than that of the biplane, which at the 
same total lift has the same induced drag as the biplane. 


APPLICATIONS OF MODERN HYDRODYNAMICS TO AERONAUTICS. 


This last remark leads us to apply also to the biplane the deduction formulas obtained for 
monoplanes. All that is necessary is to replace the biplane of span d, by a monoplane of a 
somewhat greater span kb,, which with reference to the drag—and, on the whole, with reference 
to the angle of attack—is equivalent to the biplane. If again we pass from the lift and drag 
to their coefficients c, and cy, the formula connecting the drags of any two lifting systems 
1 and 2 is 


en an Lee (60) l 
2 m Nb aba) 


in which, as is easily seen, the factor k, for a biplane having the most favorable distribution 
of lift, is the reciprocal of the square root of the second factor in formula (59). 

The tests of this formula with biplanes have shown that, when by giving a special shape to 
the wing the lift distribution was made elliptical, there was good agreement with the calcula- 


tions from monoplane experiments; with biplanes having the usual square-cornered wings,’ 


on the other hand, there was a discrepancy, which is to be attributed to the fact that the lift 
distribution on.these biplanes deviates too far from an elliptic one. We can, however, retain 
the same transformation formula if the fac- 


tor k is determined empirically for every o 7 

wing system; it is found to be somewhat 4, 

smaller than according to the theory given 

above. The experiments on this point are os a BEE 

not yet completed, so more accurate values 

can not as yet be given. The earlier Göt- 4 07 E 

tingen experiments were worked up by Dr. | | 

Munk, to whom this last idea is due, in the 05 i —— 

paper ‘‘Contribution to the aerodynamics g 5 

of the lifting parts of airplanes” in. the ae PN rl um 

Technische Berichte, Volume II, page 187. 04 _| PF IIPS |— 
27. Inthe previoussection Ihave treated 

the problem of finding the minimum of the 93 

induced drag of a multiplane, under very 9, 

definite assumptions concerning the distri- 

bution of lift over each separate wing. The 0, ai 

strict minimum problem is however differ- f PEA os 04 es 

2 


ent, viz: da 
To determine for a given front view of Fi. 54.--Curve of the function o for biplanes whose span ratio, Bun is varied. 


a lifting system that distribution oflift over 
all the lifting elements which will make the induced drag a a minimum for a specified total lift. 

In this statement of the problem the expression ‘‘a given front view”—i. e., more exactly 
stated, a given projection of the lifting system upon a plane perpendicular to the direction of 
flight—is used to mean that the wing chord is of secondary consideration, and does not need 
to be determined until later when the selection of suitable angles of attack is made. 

The general solution of this problem was also given by Dr. Munk. It will be deduced here 
in a simpler manner than that given in Munk’s dissertation, where the solution was obtained 
by the calculus of variations. By means of the stagger theorem mentioned in section 25 the 
wing system will be referred back to the corresponding nonstaggered system. For this, as we 
showed, the relation W,,= W,, holds. We shall now introduce—with the simplifying assumption 
that all the lifting elements are parallel to each’ other—a variation of the lift. distribution by 
adding at any one place a lift 6A and at the same time taking away an equal amount at some 
other place, so that on the whole the lift, which is prescribed, remains unchanged. We must 
now consider the change in the induced drag caused by this variation. If there is superimposed 
upon the lift distribution an additional air force ôA distributed. over a short portion dz, there 
arises therefrom, i in addition to the drag proper of the added lift—which, however, if sufficiently 
small is of the second order—a mutual drag, because on the one hand the added lift finds itself 
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in a flow having the downward velocity w, which is due to the lifting system, and on the other 
the lifting system is in the velocity field of the added lift. The first of these two drags, as is 


easily seen, is 5A; the other part, according to our theorem, has the same value; so the total 


drag is twice this. What condition, now, must be satisfied by the sum of the mutual drags 
caused by our twofold change of the lift distribution in order to obtain the absolute minimum 
of the induced drag? The answer is, evidently, that we will have the minimum only if by no 
change of this kind can the drag be further diminished. The sum of the induced drags, there- 
fore, can in no case be negative; also it may not be positive, because in that case by a reversal 
of the signs of the changes which we selected we could make the sum negative. Only the value 
zero is therefore allowable. Hence, if w is the vertical velocity at point 1 and w, that at point 
2, we have the relation , 

ôA, T + ALT =0 
and, therefore, since ôA, = —6A,, 

W= Wz. 

Since this holds for all the lifting elements, we have obtained the answer. The lift distribu- 
tion which in the given wing system, for a specified total lift, causes a minimum of drag is that 
which leads to the same downward velocity at all the lifting elements. With monoplanes the 
elliptical lift distribution leads to a constant downward velocity w. We recognize from this 
that the elliptical distribution in fact is that distribution of lift which causes the least drag for 
a monoplane. 

The theorem can, besides, be extended easily to the case of nonparallel lifting elements 
lying in a transverse plane. If w, is the velocity in the transverse plane perpendicular to the 
lifting element and e is the angle between the direction of w, and that of the given total lift, 
then, as may be shown without difficulty, w.=w, cose for all the elements. (If e=o, and 
hence cos e=1, the statement made above again appears.) 

28. A way to solve the problem of finding the lift distribution for a prescribed distribution 
of the vertical velocity has been indicated already in section 23. The velocity field left behind 
in the air by the lifting surface is, approximately, according to the remark made before, a 
uniplanar flow around the vortex system produced by the lifting system in its motion, and this 
vortex system may be regarded, as a first approximation, as a solid body in the fluid.. In the 
minimum case this figure, according to the results of section 27, moves like a rigid body, not 
alone in the case of parallel lifting elements, but also in the general case, for the general mini- 
mum condition, w,=w, cos e, expresses directly that the normal velocity of the fluid at an 
element of the rigid figure moving in the direction of the lift coincides with the normal com- 
ponent of the velocity wo of the rigid figure itself.. The problem is thus reduced to a perfectly 
definite one treated in the hydrodynamics of uniplanar fluid motion. 

This uniplanar flow can be brought into relation, in a specially clear manner, with the 
pressure distribution existing on the wing system. The wing system, during its motion along 
its path, imparts to one portion of the air after the other the velocities which we have learned 
to know as the result of the vortices flowing off from the wings. This transmission of velocity 
is the result of the spreading out of the pressure field of the wing system over the air particles 
one after another. In order to simplify the phenomenon for ourselves we can now imagine 
that these velocities are produced at the same moment by a sort of impulse phenomenon over 
the whole path of the lifting system. To produce this impulse it is necessary to have a solid 
figure of the shape of the geometrical region passed over by the lifting system (i. e., of the 
shape of the vortex surfaces which it leaves behind). If we are concerned with a system of 
least drag, this figure moves as a rigid body; otherwise it would also experience a change of 
shape due to the impulse. The final velocity w* of the figure coincides with the motion of the 
vortex surfaces at a great distance from the lifting system, and is therefore to be put equal 
to 2w. For a.monoplane having elliptical distribution our. figure is therefore an infinitely 
long flat plate of the breadth 6. 
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By the production of the velocity during the impulselike acceleration an increase in the 
pressure p, arises below the plate (in the case of multiplanes, under each of the plates corre- 
sponding to the separate wings) and at the same time a decrease in pressure p, above the plate 
(or the plates). We can now compare in a very simple manner the total action of the pressure 
differences at each point of a plate during the time of the impulse with the total action of the 
pressure differences of the wing in its forward movement at the point of the medium in question. 
If the resulting motion is the same in both cases, then the pressure differences integrated through 
the proper times must have the same values. If in the impulse phenomenon lasting a timera 


portion of the fluid of length 7 is considered, and if therefore the action of the lift ac in the 


time t=, required to pass over the length 7 is to be compared, the following relation must 


hold for the conditions on a strip of width dz: 
7 dA l 
ids | (p,—p,)dt= “4p de p (61) 


A formula connected with our previous relations can be obtained by a trensformation of the 
left-hand term. According to a known extension of the Bernouilli equation for accelerated 


motion we have 
of pw? 
Pit y PESO. 


For our impulse phenomenon the arbitrary time function f(t) is a constant, since at the 
points of the fluid lying far away from the impinging plate the pressure does not change. If 
the impulse is sufficiently quick, then during the short time of ee r the acceleration and 


the pressure differences will be very large, and therefore the term er may. be neglected in 


comparison with the other two, since it itself does not exceed moderäte values. We obtain 
therefore the simplified relation 


OP | p= const. = po 


p 
ot 
which, if at the beginning everything is at rest, (®,= 0) may be integrated to 
o= [m pi (62) 


we can therefore write, in equation (61), the expression p(®,—,) in place of f- Padt. The 


potential differences ,—, which here appears is, according to the connection between poten- 
tial and circulation (see sec. 5), nothing but the circulation T for a closed curve which passed 
around one edge of the vortex ribbon and intersects our vortex ribbon at z, the point consid- 
ered. This circulation is again nothing but the circulation around the wing at the point z. 
Tf the factor ldx is omitted from both sides of equation (61), it takes, as a result of this trans- 


formation, the form 


dA 
dt = p(®, — 8) V=pVI (63) 


We have thus proved in an entirely independent way, as we see, the Kutta-Joukowski 
theorem for a wing element, which previously we took over, without proof, from the infinitely 
long wing. 

The relations deduced in the previous paragraphs permit, in the case of a constant w, 
the formation of general theorems for w and Win place of (39) and (40). By integration 
of (63) the total lift is at once obtained 

=pV2f(®,—%,)de (64) 


The values of ® in this formula are proportional to the velocity w* of the vortex ribbon, 
that is, are dependent upon A. Quantities which are independent of A are derived if the poten- 
tials are divided by w*. In this way we obtain the potentials for a velocity of the vortex 
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ribbon equal to 1. The potential, being the line integral of the velocity, has the dimension 
velocity times length; the potential & for w*=1 has therefore the dimension of a length, and 
hence 


Zf% -H)de 


is a surface, which will in what foilows be called F’, which depends only upon the geometrical 
properties of the projection of the wing system upon a plane perpendicular to the direction 
of flight, therefore upon the front view of the wing system; and which evidently for geometrically 
similar front views is proportional to the square of the span. By introducing F” into equation 
(64) we have, since = ¢u*, 
A=pVw* F" (65) 
From this we may immediately deduce w*, and thereby also the downward velocity at the point 
of the wing system 
En A 
w= 5w “Dp VF’ (66) 


If this value is introduced into the relation W= PA, we have 
A? A? 
= op V2F" yQ (67) 
The evaluation in the manner of the flow of figure 49 gives a potential 4, if the span of 
the wing is set equal to b, which has the value (6/2)?—2? at the plate. The geometrical expres- 


W 


: 2 
sion of this value gives a circle having the span b as-diameter, therefore pate. Using this 


value formula (67) passes over in fact into formula (40). 

It may also be noted that a uniform velocity can be superimposed upon the uniplanar 
flow here discussed, whose discontinuity in potential at the rigid figure representing the vortex 
ribbon causes the surface F’, without thereby changing the relation for F”, for the potential 
discontinuity between the lower and upper sides, with which we are here concerned, is not 
changed by the superimposed uniform motion. We may now choose the velocity of the uni- 
form motion exactly opposite and equal to the velocity w* of the rigid figure, and thereby 
secure the condition that in the new flow the rigid figure is at rest and is surrounded by a flow 
which .at.infinity has the velocity w. The forces which the rigid figure experiences by the 
production of this motion, and which are connected intimately with the so-called “apparent 
mass,” are what. we have here set in parallel with the wing-lift. 

. The surface F’ supplies in addition a very simple mechanical connection between the 
velocity w on the one hand and the lift and drag on the other. According to equation (65) 


A=pF’ Vw* 
2: 
WV = Aw=pF'V 


where in the second equation use has again been made of the relation w*=2w. Now pF” V 
is the mass of air flowing per second through the section Æ”. If in order to simplify the whole 
problem it is once assumed that all the air particles within the section F” experience the full 
deviation w* but that all outside are entirely undeviated, then exactly the correct lift and 
the correct work due to drag are obtained by application of the impulse theorem and the energy 
theorem. For the lift.is now equal to the mass of the fluid deviated per second times the 
vertical velocity imparted to it, therefore equal to the impulse imparted to the medium. Also, 
in the same manner, the work done per second by the drag WY is the product of the mass of 
the fluid passing per second times the half square of the deflection velocity, and therefore equal 
to the kinetic energy left behind in the medium. This relation is indeed best suited to establish 
the phenomena of the theory of airplanes in a course for students who are only slightly skilled 
in mathematics. The fact that for a monoplane the circle having a diameter equal to the span 
comes out as the surface F” is one that will appear most plausible to the laity. 
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29. Of the theory pictured in the preceding section, according to which the determination 
of the induced drag in the case of the most favorable lift distribution is reduced to a problem 
of the potential theory, manifold applications have already been made. Especially, Dr. Gram- 
mel and K. Pohlhausen have treated, at my instigation, the case of the biplane made up of two 
straight monoplanes of the same span, and also, on the other hand, that of a monoplane having a 
longitudinal slot. The calculations in both cases are solved by means of elliptic integrals. I 
have given the formulas in my Wing Theory II. It may be sufficient here to state the practical 
final result, which is referred to the magnitude of the surfaces F’. These surfaces are best 
expressed for biplanes in terms of the corresponding surfaces of the monoplane having the same 


span. In fact, the ratio F” : A b2, as is easily seen, equals the square of the factor k, introduced 


in section 26, by which the span must be increased in order to have a monoplane of the same 
induced drag. The values of k? for the biplane are obtained from the following table. The gap 
of the biplane, i. e., the distance apart of the two wings, is designated by h. 


TABLE 2, 


Values of k? = F’ if pe 


0,05 


. 0.10 
1.156 


1.212 


0,15 


0.2 
1.289 1.3 


0 | 0.05 | 010 | 015 | 02 103 | 0A fasan 
1.000 


hjb=... 0.3 
| Rs 52 1.461 


The values given in the table may be expressed by the approximation formula 


1.027 +3.84h/b 
TT 1.68h]6 (68) 
In the case of the monoplane having a slot a suitable comparison wing is obtained by shov- 
ing the two halves of the monoplane together until the slot is closed. If bis the original span 
and d is the width of the slot, this monoplane has evidently the span d—d. We shall therefore 


T 


form the ratio Æ”: 7 (b—d)? and again designate it by k. Calculations gave the following 
values: 


k 


TABLE 3. 


Values of k? = Fr: (bd)? 


| 
0. 0316 | 0.100 


0.000 
0. 6200 | 0.568 


-| 1.000 


0.001 
0. 762 


0.010 
0.676 


0.250 
0.528 


0.500 
0.506 


1.000 
0.500 


It is seen that even very narrow slots produce an important increase in the induced drag. 

For a very wide opening i? falls to one-half, as may be deduced easily from the fact that now, 

instead of one monoplane, we have two monoplanes of half the span. The values given in the 
. table may be expressed by the approximation formula 


1 
 2y1+0.35 (log,, b/d)? (69) 


Figures. 55-57 show, at the left, the uniplanar w* — flow and, at the right, the surfaces F”, for a 
monoplane, a biplane, and a monoplane with a slot. 


=I 


F. AEROFOILS IN A TUBE OR IN A FREE JET. 


30. To draw conclusions from the experimental results obtained in a tube bounded by 


solid walls or in a free jet from a nozzle, it is very useful to know the influence of the neighbor- ` 


ing walls and of the boundaries of the jet upon the phenomena at the aerofoil. We wish indeed 
to know the behavior of the aerofoil in an air space infinitely extended in all directions; and the 
problem therefore arises to introduce a method for passing by calculation from the case which 
prevails in the experiments to that of the unlimited air space. For this purpose we shall next 
state clearly the boundary conditions which exist at solid walls parallel to the direction of the 
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FIG. 57, 


Left: Flow behind a monoplane, a biplane, and a monoplane with a slot, with reference to an 
observer moving downwards with the velocity W* of the vortex system. 
Right: Corresponding surfaces Fi, 
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wind and at the free boundary of a jet. At solid walls, the velocity components normal to the 
wall,w,, must equal zero; on a free jet boundary, on the other hand, the pressure is to be put 
equal to that of the surrounding quiescent air layer, and therefore is constant. We can trans- 
form these last relations as follows, keeping within our theory of the first order. According to 
the Bernouilli theorem, if V is the undisturbed wind velocity, and u, v, and w are the additional 
velocities 


p+ S| w+ (rosa |=p+ $V 
or, since p= Po, 
W+e+w+2Vv=0. 

If we neglect the squares of the disturbance velocities as being small of the second order, we 
have as the approximate initial condition for the free jet v=o. We proceed a step farther upon 
the path indicated to us by the approximation theory of the first order if we prescribe the value 
v=o, not for the actual jet boundary, but for that cylinder which is given by the’ surface of the 
undeviated jet. By doing this the boundary condition for the free jet becomes very similar, 
in a formal way, to that for the solid walls. 

The two problems can now be solved in the following manner: We consider first the velocity 
field for the unlimited air space, according to the exposition previously given. This field offers, 
both in the case of the tube and in that of the free jet, contradictions with our boundary condi- 
tions at the walls or the jet boundary. We must superimpose a velocity field which in the 
interior of the region considered is free of singularities and which on the boundaries has veloci- 
ties opposite to those velocity components, the vanishing of which is prescribed by the boundary 
condition. It is easily seen that by the superposition of this second velocity field on the original 
one the boundary conditions are satisfied exactly. The influence of this second field upon 
the aerofoil is now exactly that influence which we are seeking, and which we can calculate 
from the results of the theory of aerofoils as soon as this second field is known. 

The additional velocity field corresponds to a pure potential motion; we have, therefore, the 
problem of determining its potential 6. In the case of solid walls we are thus led to the problem 
of finding the potential for a given region (the interior of the tube) when the normal component 
Wn of the flow is given at the boundary of the region. This is the so-called ‘‘second boundary 
value problem” of the potential theory. The corresponding problem for a jet, as we shall see 
at once, leads to the ‘first boundary value problem,” in which at the boundary the values of 
the potential itself are prescribed. According to what has been said above our region is a cylinder 
whose generating lines are parallel to the velocity V, hence, parallel to the Y axis, and for each 


point on the boundary the relation 3 — (in which the dashes indicate boundary values) is 


prescribed. Integrating this relation for each generating line gives 


rd y - 
& (y= he ody 

If we go sufficiently far upstream every influence of the aerofoil vanishes; therefore for y= —- œ, 

=o; and hence y= — œ is taken as the lower limit of the integral. By this, then, we obtain 

the boundary values of the potential ® (y). 

The complete calculation of the added potential ® for the entire interior of the tube or jet is 
fairly difficult. If we concern ourselves, however, only with our main problem,. to determine 
the corrections which must be applied to our experimental results, then we can again assume 
that the velocity components perpendicular to the axis of the tube in the plane of our aerofoil 
are half as large as at a great distance behind it. This consideration, which proceeded from the 
comparison of a vortex filament proceeding to infinity in one direction only with one proceeding 
to infinity in both directions, holds here exactly as well as in the cases discussed previously. 
We can therefore pass here as before from the space problem to a uniplanar one if we calculate 
the phenomena far behind the aerofoil. Our boundary conditions for the uniplanar problem are, 
for the tube, Wn =0, for. the free jet, ®=const. The last condition may be interpreted specially 
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conveniently if use is made of the method of treatment of section 28. Since in this, for the end 
of the impulse phenomenon, = f (po— p) dt, ®=const. means simply that p= const., which, 
indeed, was the original boundary condition for the free jet. 


31. The conditions stated in the preceding section can be secured most easily for a jet, or 
tube, of a circular cross section. In this case the added motion is obtained very simply by 
assuming for every vortex flowing off an equally strong one outside the circle, at the point out- 
side corresponding to the one inside according to the reciprocal radii. If the direction of rota- 
tion of the external vortex is taken the same as that of the interior one, then at the points of 
the circle the boundary condition for a free jet is obtained; and, if opposite directions of rotations 
are taken, then the boundary condition for a tube is satisfied. This may be expressed by saying 
that there is combined with the aerofoil another obtained by reflexion according to reciprocal 
radii, whose circulation at corresponding points is the same in absolute value as that of the 
actual aerofoil, and for the jet it has the same sign, but for the tube the opposite sign. 

The exact calculation has been made for a straight monoplane in the middle of the jet, 
assuming the lift to be distributed according to a half ellipse. If b is the span of the monoplane 
and D the diameter of the'jet, then the disturbance velocity w’ caused by the jet boundary at 
the distance x from the middle of the jet is 


,. A 3,,5 35 
w vl! +IP + get a tete.) (70) 
in which ¿= 2gb/D?. 
The added drag calculated from this velocity according to equation (38) is found to be 


‚___ 4 3(b\, 5/5. 
w - m +) gD) oe | ae 


A similar. calculation for a uniform lift distribution gave for the first term in the formula 
for the drag the same value as in equation (71). It appears that the other terms of the series 
have but little importance with the usual ratios, so that we can limit ourselves to the first term. 
An approximation treatment shows, further, that any small wing system, in the middle of the 
circular jet gives rise to the same expression. We can therefore write for the total induced 
drag of the wing system in a jet of cross section Fù, if the surface F” is again introduced from 


section 28, 
42/1 1 
ar Ge oF A) : (72) 


For a tube of circular cross section the same disturbance effect is found, but with the oppo- 
site. sign; and therefore we have the approximation formula for the drag 


A?/ 1 1 
ae ar) (72a) 
The correction, owing to the consideration of the finite cross section of the jet, is for the 
ratios ordinarily used not small. For 2-4, it is already one-eighth of the induced drag. 


Formula (71) gives 0.1262 instead of 0.125; the corresponding formula for uniform distribution 
gives 0.127. Itis seen, therefore, that the differences are not great, and that the approxima- 
tion formula (72) is satisfactory for most cases. 

For a tube of rectangular cross section the calculations would have to be made in such a 
manner that the aerofoil was mirrored at all the walls an infinite number of times, like a check- ` 
erboard. Further development of the calculation leads to elliptic functions. It has not yet 
been carried through. One can assume, however, that for a tube having a square cross section 
the influence of the walls will be of about the same magnitude as for the circle having an equal 
area. : 


APPLICATIONS OF MODERN HYDRODYNAMICS TO AERONAUTICS. 


G. APPLICATION OF THE THEORY OF AEROFOILS TO THE SCREW PROPELLER. 


32. The fundamental ideas of the aerofoil theory can be applied step by step to the screw 
propeller. For the elements of the blades the Kutta-Joukowski formula holds, viz, that the air 
force is perpendicular to the velocity c of the element with reference to the air and-that, per 
unit length of the blade, it has the value pre. Corresponding to what has gone before, vorti- 


ces will arise at the blade, having a vortex strength per unit length equal to = - If we wish 


again to construct a theory of the first order, that is, if we agree to consider as small the air 
forces and the velocities produced by them, then again the proper motion of the vortices will be 
small and therefore in a first approximation may again be neglected. The vortices then have 
the shape of screw lines and form vortex ribbons which—if for the sake of simplicity we assume 
straight radial blades—have the shape of ordinary screw surfaces. 

- The calculation of the velocity field of a screw vortex is markedly more complicated than 
that of a rectilinear vortex and leads to functions which thus far have not been studied in detail. 
In spite of this it is possible, as Dr. Betz has shown, to prove a series of general theorems very 
similar to those of Munk for multiplanes. Since the velocity ¢ is not the same at the separate 
blade elements, we must speak of the “ work lost” where Munk speaks of drag. The work applied 
for the motion of the propeller is composed of two parts—useful work + work lost. The latter 
in our ideal case, where friction is excluded, is transformed completely into kinetic energy of the 
air. The kinetic energy stands again in close connection with the vortex system produced by 
the propeller. Betz proved, among others, the following theorems: 

(1) If two elements of a propeller blade lie upon the same radius at distances x and ¢ from 
the axis, then the work lost at the point & due to the disturbance velocity caused by the air 
force at the point x is equal to the work lost at the point x owing to the disturbance velocity 
caused by the air force at the point £. 

(2) This theorem must be somewhat modified for two elements which do not lie on the same 
radius. It reads: The work lost at the point & due to the disturbance velocity caused by the 
air force at the point x is of the same amount as the work which would be lost at the point x if 
the screw vortex proceeding from the element at £ were to pass out forward in the prolongation 
of the actual vortex instead of going backward.” 

(3) This last theorem leads at once to the following relation for the sum of the two amounts 
of work lost: The total work lost due to the mutual action of the air forces by two blade elements 
at points x and é is the same as the work which would be lost at one point alone if the screw vor- 
tex proceeding from the other point were to extend to infinity both forward and backward. 

It is easily seen that this theorem is perfectly analogous to the stagger theorem of section 
25, for if the vortex of the inducing element extends in both directions, then the position of the 
element itself on its own vortex strip is immaterial as far as the velocity field produced is con- 
cerned. It is therefore true of screws that nothing is changed in the total energy-loss if blade 
elements are displaced in any way, without change of their air forces, along the relative stream- 
lines passing through them (i. e., in this case, screw lines). This naturally is connected again 
with the fact that the total amount of the energy loss depends only upon the final distribution 
of the vortex systems, not upon the relative position of the places where the separate vortices 
arise. 

Theorem No. 3 will be of use to us also in what follows. It can be made clearer by the fol- 
lowing consideration. The field of the vortex ribbon of a lifting element dies away very quickly 
forward of the element, but in the rear it extends over the entire length of the path traversed. 
If the sum is formed of the two mutual losses in work of two elements at the points x and £,we 
can proceed, owing to the stagger theorem, to displace one of the two elements along its screw 
line so far backward that its velocity field is no longer appreciable at the position of the other 


12 In this the sense of rotation of the transverse vortex is to be reversed. 
13 The transverse vortex in this case cancels out completely in the determination of the velocity field, since it appears twice with opposite 


senses of circulation, 
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undisplaced element. At the same time, however, the influence of the latter element upon the 
first is increased since its vortex ribbon, viewed from.the new position of the first element, extends 
as far forward as backward. The sum of the two mutual losses in work is reduced in this manner 
to the loss which the velocity field of the front element produces upon the displaced one." 

(4) The most important of Betz’s theorems, from a practical standpoint, furnishes the com- 
plete analogy to Munk’s theorem concerning the wing system having. the least drag, and, corre- 
sponding perfectly to the statements in sections 27 and 28, may be expressed thus: The flow 
behind a propeller having the least loss in energy is as if the screw surfaces passed over by the 
propeller blades were solidified into a solid figure and this were displaced backward in the 

nonviscous fluid with a given small velocity. The potential differ- 
ence between the front and rear sides of a screw surface at one and 
the same point furnishes, then, again the circulation T of the corre- 
sponding point of the propeller blade. 
A short proof of theorem 4 will be given. For this purpose the 
Fo principal equations for the action of a screw must first be deduced. 
The screw is imagined to be displaced with the velocity v relative 
to the air, and to rotate at the same time with the angular velocity w. 
A blade element at the distance x from the axis has then, with ref- 
erence to the air which in the theory of the first order may be 
assumed to be at rest, the velocity c, with the components v and aw. (See fig. 58.) 

If no vortex were produced, then, with the assumption of a nonviscous fluid, an air force 
dP would arise, which, according to the Kutta-Joukowski theorems, would be perpendicular 
to the velocity ¢ and would have the value, for a blade element of length dz, 


dP 


da 
w 


Fig. 5S.—Velocity field near a blade 
of a screw-propeller. . 


dP = pV cdz (73) 


The force dP is decomposed into two components, of which the one in the direction of v 
interests us specially, since it is applied to the screw. This component is 


dS=dP . cos e=pTxwdz (74) 
The total thrust, if there are n blades, is then 
S= por |T x de (75) 
The other component 
dT=dP sin e=pTv dz (76) 


furnishes a contribution as a torque to the rotation moment. It is seen at once that dS . v= 
dT . cw, i. e the useful thrustwork is equal to the work done by the torque hitherto used 
in our calculations. This depends immediately upon our assumption that the force dP is 
perpendicular to the velocity c. But the screw blades actually produce a vortex system and 
we must ask as to the reaction of the vortex system upon the phenomena of a screw. We 
shall assume, exactly as in the aerofoil theory, that we turn the blade profile in such a manner 
that the lifting forces desired by us actually come into play. Since we are interested here 
merely in the loss in work caused by the vortex system, we have to do only with the drag 
components caused by the vortex system. This depends, exactly as before, upon the velocity 
component perpendicular to the velocity of motion of the element, which in this case equals c. 
We shall again designate it by w. The added velocity component w furnishes a drag in the 
direction of motion equal to 


dQ= dP “= lw dr (77) 
The loss of work per second is therefore 
dQ .c=p¥de . w . c (=dP . w) 


4 This process of thought can be applied, naturally, in the same way to aerofoils and furnishes a convenient deduction for the sum of the 
drags Wirt Wa. 


APPLICATIONS OF MODERN HYDRODYNAMICS TO AERONAUTICS. 


If now, according to figure 58, we put c=v/sin E, our problem is to make a minimum the total 


loss of work 
n r T 
L=pvS Í na (78) 
1 


The variation of this quantity must therefore be put equal to zero. We shall proceed, for 
this purpose, similarly to the way Munk’s theorem was deduced in section 27. We shal 
change by small amounts the circulation at two places, which may reach from x, to #,+ da, 
and from z, to &,+dx,, in such a manner that the total thrust remains unchanged. According 
to equation (75) we must make 

ôT, . xde, +ôr, . zdz, =0 (79) 
Exactly as before the condition for the minimum is obtained if the loss of energy due to our 
added circulation remains unchanged. In order to calculate the loss, let us make use of 
theorem No. 3 and assume that the added wing forces are brought into action.far behind the 
propeller so that the loss is merely the product of the added air force by the velocity w* 
which arises from the vortices of the propeller, and therefore for the first element’ is equal to 


* eat 
pvöT, Ih dz, Omitting the constant factor, we obtain as the minimum condition— 


w,* 
dtl) 3 = ne =: 
from which is derived, making use of pantie (79) 
* DR 
u inet: (80) 


Tı sine 2%, SiN e, 

We must compare this condition with that obtained for the velocity components normal 

to a rigid screw surface, when this surface is moved backward with the a 
velocity w’. We then have (see fig. 59): i 


Wn =W COS €. 


But on a screw surface the pitch A is connected with the angle of piteh € 
and the radius x by the relation 4=27z tan e. 
Multiplying this last equation with that for wa, and pone we get- 


2rw! £ sin € 
Un a 


and therefore 


= _ _ const. (81) 


On comparing (81) with (80) it is seen that by a suitable choice of w’ the 
value of wa can always be made to agree with that of w*, which proves F14. 50—Ideslived vortex sys- 

; tem of a screw-propeller. 
Betz’ theorem. 

33. In order to learn more accurately the nature of the distribution of circulation which 
we are seeking, we shall proceed as if the velocity field at great distances from the screw is pro- 
duced by having the velocity w’ in the direction of the axis imparted impulsively to the rigid 
figure composed of the screw surfaces. In a purely qualitative way one can see that with any 
system of screw surfaces having a small pitch the air in the interior of the system is actually 
accelerated backward, with, of course, the appearance of tangential velocity components whose 
intensity is a function of the angle e and is greatest fore=45°. At the axis itself there is neither 
an axial nor a tangential necaloration:. ines simple are the conditions near the outer boundary 
of the screw surface where a flow around the edges of the surfaces occurs. 

In order to obtain a quantitative statement, we shall for simplicity’s sake next think of a 
screw having a large number of blades. Our rigid figure consists, then, of a very large number 
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of screw surfaces, lying close together, and therefore the air is led with difficulty into the interior. 
When the impulse occurs it escapes in the direction of the normals to the screw surface. The 
radial velocity components w, will be appreciable only in the neighborhood of the outer boun- 
dary of the screw surfaces; further in, we may put it equal to zero approximately. For the 
tangential components w, and the axial components Wa, the relations hold, as is easily seen, 


w,=w* sin e=w’ cos esin e 
Wa = w* cos e=w’ cos? €. 
The angle may. be expressed by writing 
r 
tan = 5 ==, (82) 


In this, for brevity’s sake, hat is put equal to r’ (r’ is that radius for which the pitch of the 
screw, tan e =1). Then 


i 


” er and cos e= iF Ta (83) 
and hence 
; re z 
wmi agg ANE nw Ze (84) 


We must now determine the circulation around the separate blades as a function of the 
radius x. For a screw with n blades the total circulation of the vortices inside the circle of 
radius x coincides with the line integral for the closed circle of radius x; this circulation must 

evidently equal nT, where T is the circulation of one of the screw blades at 
| the point x. From this we have 
In» Wy _ 2rr'w' x? hs 
aa ae ee aa > 
The curve for T according to equation (85) is shown in curve I of figure 62. 

At the ends of the blades we would expect to have a decrease of the cir- 
culation of the same character as found for aerofoils. An approximate treat- 
ment can be devised in the following way: We imagine an infinite series of 
aerofoils which have a distance apart a and are not staggered and which 
extend infinitely far toward the left. We inquire what is. the most favor- 
able distribution of lift near the ends of these aerofoils. The distance a is 
| then to be made equal to the perpendicular distance apart of the edges of 
two consecutive blades of the screw, i. e., according to figure 59, 


Fis. 60.-—Conformal trans- 


formation , 


@ 1/1 _h _ 2ar r? 
z= logs (13). am, en zn F7 


(86) 
The problem now, according to the procedure of section 28, may be solved by seeking the po- 
tential flow around the edges of the corresponding family of planes and by determining the dis- 
continuity of potential at the planes. This problem may be solved without difficulty by means ~ 
of conformal representations. (See sec. 10.) It can be shown that the plane with the straight 
cuts as shown in figure 60, which we shall call the z plane, may be transformed into the unit 


circle (£ plane) by the formula 
758 t 
The flow of figure 58 is transformed thereby into circulation flow around the unit circle; in fact 
$44 V=10 logt 
After a short calculation, by elimination of t, we have 


a +1 
z= log cos ~% (88) 


APPLICATIONS OF MODERN HYDRODYNAMICS TO AERONAUTICS. 


For the surface around which the flow takes place, and which is given by the streamline y = 
we have, therefore, 
z= log cos 4 
ie C 
or 
= + C cos ea (89) 


which gives the real values for negative values of x. The potentials thus obtained or the 


velocity equal to 1 of the free flow (to obtain which O must be put equal to 2), which, accord- 


ing to what has gone before, give us the surface F”, form a picture such as is shown in figure 
61. By means of this one can form a definite judgment as to how the circulation, and with it 
the thrust also, decreases at the blade tips. We can replace the shaded portions of Figure 61 
by a straight line, having an equal area below it which, in accordance with , 

the integration performed, must lie behind the blade tips at the distance 


a' =a log g Ž=0.2207a (90) mm 
We conclude from this that, with screws also, the decrease of circulation 
at the blade tips has about the same effect as if the screw had a radius dimin- 
ished by 0.2207a, and then the air would be considered uniform in every circle of I 
radius x (as would be the case for a screw having an infinite number of blades). ra. 81. Potential ob- 
The properties found for. the inner portion of the screw and for its edge tained by the flow of 
may be combined into a single formula, which can be applied as an approxima- *#""°® 


tion formula also for screws having a small number of blades. This formula is obtained by 
== (fx) 
multiplying the value in equation (85) by the expression — = cose “a” , which for large values 


of” = takes the value 1. (For —z of formula (89) r—z is here substituted, as is obvious.) 


Thus we obtain the formula 


Pah 2 —r (T—X) 
Sa Eg ote (91) 
The curve of T according to equation (91) is given in figure 62, fora 4-blade screw and for 
r’:r=1:5, which correspond to average conditions in practice. 
The whole deduction holds, as has already been remarked, for screws which are not heavily 
loaded. For screws with heavy loads an improvement can be intro- 
duced by calculating the pitch of the screw surfaces formed by the 


T= 


Asymptote _A_ 


r 
j vortices, corresponding to the state of flow prevailing in the circular 
C 3 [$ iti ar i 
xr’ ic. plane of thescrew. Instead of writing tan = zo’ we must write, more 
a 
y+ ve 
Fig. 62.—Distribution of circu- : y A i : : è 2 
lation, I. Infinite number of exactly, tan & = ‚in which V is the velocity of flight, since in 
blades. 1I. Four blades « is Lu — at 
the distance of two vortices. 2 
eens the screw disk plane half of the final disturbance velocities is already 


+ 
present. A useful approximation is obtained if, retaining our formulae, v is put equal to V+5 i 
w’ 


and therefore r’ is put equal to V+ 2. 
w 


After the circulation is known, the distribution of thrust and torque may be calculated 
easily by means of equations (74) and (76), and thus, following the method used in the aerofoil 
theory, the requisite widths of the blades and angles of attack may be determined in order that 
for a given working condition (i. e., r’ and w’ given), in which the screw is to have the most 
favorable performance, all the information may be deduced from the theory. By taking into 
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account the more exact velocity relations in the propeller-disk plane this information may 
be improved. 


The aerofoil theory has numerous further applications. An investigation of curved flights 
specially of themoment—important in discussions of stability—around the longitudinal axis in the 
case of a wing moved in a circle, is at present being made, also the calculation of the moment 
of a warped wing. A series of not unimportant single questions must wait for a further im- 
provement of the theory, e. g., various conclusions specially concerning properties of profiles, 
influence of curvature, etc., can be reached, if we pass from the lifting line to the case of a load 
distributed also along the chord; for the treatment of a wing set oblique to the direction of 
flight the assumption of a load distributed along the chord is necessary since in this case the 
conditions contradict the “lifting line.” ‘Investigations of this kind, which can be accom- 
plished ‘only by very comprehensive numerical calculations, were bean during the war but 
since then, owing to a lack of fellow workers, have had to remain unfinished. A similar state- 
ment also applies to'the calculations of a flapping wing already begun, in which one is likewise 
forced to assume the lift distributed along the chord, since otherwise the result is indefinite. 
Therefore much remains to be done. 


MOST IMPORTANT SYMBOLS. 
p =density. 
14 =velocity of the airplane. 
u, v, w=velocity components in the X, Y, Z directions. (In the case of an airplane X is in the direction. of the 
span of the wings, Yis in the direction of flight, Z is vertical.) 


q= eV =dynamical pressure. 


b =span of a wing (“Breite”). 

t =chord of a wing (“Tiefe”). 

h =gap of a biplane (“Höhe”). 

£ =area of surface (=b . t) (“Fläche”). 
A =lift (“ Auftrieb”). 

Ww =drag (“Widerstand”). 


= fmt coefficient (=2 K, “‘absolute”). 


= A =drag coefficient (=2 K, “absolute ”). 


a =angle of attack. 

T =circulation. 

® =velocity potential. 
yp =stream function. 


LIST OF THE MOST IMPORTANT LITERATURE. 


Abbreviations: Z. F. M.=Zeitschrift für Flugtechnik und Motorluftschiffahrt. 
TB =Technische Berichte der Flugzeugmeisterei. 


A. GÖTTINGEN PUBLICATIONS. 


L. PranprL: Tragflügeltheorie I. und II. Mitteilung. Nachr. von der Kgl. Gesellschaft der Wissenschaften. 
Math.-phys. Klasse 1918, S. 451, u. 1919, 5. 107. (Aerofoil Theory, I and II Communications. Nachr. Kgl. 
Gesellschaft der Wissenschaft. Math.-phys. Classe, Göttingen, 1918, p. 451, and 1919, p. 107.) 

—— Tragflächen-Auftrieb und -Widerstand in der Theorie. Jahrb. der Wissenschaftlichen Gesellschaft f. Luftfahrt, 
V. 1920, 8. 37. (Aerofoil Lift and Drag in Theory. Jahrb. d. Wissens. Gesells. f. Luftfahrt, V. 1920, p. 37.) 
—— Der induzierte Widerstand von Mehrdeckern. TB Bd. III,.S. 309. (The Induced Drag of Multiplanes. TB. 

Vol. III, p. 309.) 

G. Fuurmann: Theoretische und experimentelle Untersuchungen an Ballonmodellen. Jahrbuch 1911/12 der Motor- 
luftschiff-Studiengesellschaft, 8. 65. (Theoretical and Experimental Investigations of Models of Airship Bodies. 
Jahrb. 1911/12 d. Motor Luftschiff-Studiengesells, p. 65.) 


APPLICATIONS OF MODERN HYDRODYNAMICS TO AERONAUTICS. 


A. Betz: Die gegenseitige Beeinflussung zweier Tragflächen. Z. F. M., 1914, S, 253. (The Mutual Influence of 
Two Aerofoils. Z. F. M., 1914, p. 253.) 


—— Untersuchung einer Joukowski’schen Tragfläche. Z. F.M., 1915, 8.173. (Investigation of a Joukowski Aero- 


foil. Z.F. M., 1915, p. 173.) j 

—— Einfluss der Spannweite und Flächenbelastung auf die Luftkräfte von Tragflächen.. TB, Bd.1,8.98. (Influ- 
ence of Span and Wing-Loading upon the Air forces of Aerofoils. TB, Vol. I, p. 98.) 

—— Beiträge zur Tragflügeltheorie mit besonderer Berücksichtigung des einfachen rechteckigen Flügels. Disserta- 
tion, Göttingen, 1919. Auszug in Beiheft II der Z. F. M., 1920, S. 1. (Contributions to the Theory of Aerofoils 
with Special Reference to the Simple Rectangular Wing. Dissertation, Göttingen, 1919. Extract in Beiheft. 
Il Z. F. M., 1920, p. 1.) 3 

—— Schraubenpropeller mit geringstem Energieverlust, mit einem Zusatz von L. Prandtl, Nachr. v. d. Kgl. Gesell- 
schaftder Wissenschaften, Math.-phys. Klasse 1919, S. 193. .(The Screw Propeller having the Least Loss of Energy, 
with an Appendix by L. Prandtl. Nachr. Kgl. Gesellschaft der Wissenschaft, Math.-phys. Class, Göttingen, 
1919, p. 193.) ' 

—— Eine Erweiterung der Schraubenstrahltheorie. Z. F. M., 1920, S. 105. (An Extension of the Theory of Screw 
Jets. Z. F.M., 1920, p. 105.) : 

M. Munk: Beitrag zur Aerodynamik der Flugzeugtragorgane. TB, Bd. II, 8.187. (Contribution to the Aerodynamics 
of Lifting Airplane Members. TB., Vol. II, p. 187} i 

—— Isoperimetrische Aufgaben aus der Theorie des Fluges. Dissertation, Göttingen, 1919. (Isoperimetric problems 
from the Theory of Flight. Dissertation, Göttingen, 1919.) 

C. WIESELSBERGER: Beitrag zur Erklärung des Winkelfluges einiger Zugvögel. Z. F.M., 1914, S. 225. (Contribution 
to the Explanation of the Formation Flight of Migratory Birds. Z. F. M., 1914, p. 225.) 

—— Experimentelle Prüfung der Umrechnungsformeln. In “Ergebnisse der Aerodynamischen Versuchsanstalt 1, 
Lieferung,” München, 1921, 8. 50. (Experimental Examinatjon of the Formulae for Relations between Aero- 
foils.. In ‘Ergebnisse der Aerodynamischen Versuchsanstalt,” Part I, Munich, 1921, p. 50.) 


B. WORKS ON THE TWO DIMENSIONAL PROBLEM.- 


W.M. Kurra: “Auftriebskräfteinströmenden Flüssigkeiten.” Illustr. aeronaut. Mitteilungen, 1902, 8.133. Ausführ- 
lichere Abbandlungen in den Sitzungsber. der Bayr. Akad. d. Wiss., Math.-Phys. Klasse 1910, 2. Abh. und 1911, 
8.65. (Forces of Lift in Flowing Fluids. Illustr. aeronaut. Mitteilungen, 1902, p. 133. More detailed papers 
in the reports of the proceedings of the Bavarian Acad. of Sciences, Math.-Phys. Class 1910, 2d report,. and 
. 1911, p. 65.) l l 

N. Jouxowsxı: Ueber die Konturen der Tragflächen der Drachenflieger. Z. F.M., 1910, 8.281. ‘‘Aerodynamique” 
aus dem Russischen übersetzt von Drzewiecki, Paris, 1916. (On the Contoursof Airplane Wings. 'Z.F.M.,1910, 
p. 281. “Aerodynamics” translated from the Russian by Drzewiecki, Paris, 1916. 

R. Gramuer: “Die hydrodynamischen Grundlagen des Fluges.’? Braunschweig, 1917. (The Hydrodynamic Princi- 
ples of Flight. Brunswick, 1917.) j 

R. v. Mises: Zur Theorie des Tragflächenauftriebs. Z.F.M., 1917, 8.157, 1920, S. 68 und 87. (The Theory of Lift 
of Wings. Z. F. M., 1917, p. 157, 1920, pp. 68 and 87.) 


After this memoir was written, two papers, by R. Fuchs and E. Trefftz, on the theory of aerofoils appeared, both 
of which discuss the theory of a monoplane and that of airplanes of least drag. . These papers are published in the 
“Zeitschrift für Angewandte Mathematik und Mechanik,” 1921, Heft 2 u. 3, Berlin. 


55. 


